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|_i \ A fully consistent linear perturbation theory for cosmology is derived in the pres- 

ence of quantum corrections as they are suggested by properties of inverse volume 
operators in loop quantum gravity. The underlying constraints present a consistent 
deformation of the classical system, which shows that the discreteness in loop quan- 
Q\ ' turn gravity can be implemented in effective equations without spoiling space-time 

covariance. Nevertheless, non-trivial quantum corrections do arise in the constraint 
algebra. Since correction terms must appear in tightly controlled forms to avoid 
anomalies, detailed insights for the correct implementation of constraint operators 
can be gained. The procedures of this article thus provide a clear link between 
fundamental quantum gravity and phenomenology. 

o 



1 Introduction 

Quantum gravity is expected to play a role in the early universe, in such a way that it 
may become subject to observational tests in cosmology. As is well known, a complete 
theory of quantum gravity is difficult to construct, and even if one would have a fully 
convincing candidate it would remain difficult to link such a fundamental formulation to 
clear-cut observational consequences. Daunting as this may seem, such a problem is not 
specific to quantum gravity and has been circumvented highly successfully in other areas. 
For instance, to date there is no complete and fully rigorous construction of interacting 
quantum field theories on flat space-time, and yet clear and experimentally well-tested 
procedures to extract predictions have been used for decades. 

Quantum gravity certainly does have additional problems which do not arise in quantum 
field theories on flat space-time. Paramount among these issues, related to the general 
covariance of the theory and thus to consequences of the fully constrained nature, are the 
problem of time, the self-interacting nature of gravity and the notion of a physical Hilbert 
space. Yet, quantum field theory is a lesson that the lack of a completely formulated 
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underlying theory should not prevent one from making trustworthy statements valid at, 
e.g., low energies. Not all the mathematical constructions, whose well-defined existence 
one may wish to prove rigorously, are required for such purposes. They are surely necessary 
at a fundamental level, and they have often stimulated much further research. But they 
do not directly relate to observables, and thus can, for some purposes, be ignored. 

The key tool for extracting potentially observable consequences without being para- 
lyzed by open issues in a fundamental framework are effective formulations. They capture 
quantum effects by describing relevant aspects of an evolving wave function. They allow 
one to focus on the relevant degrees of freedom, such as expectation values or fluctuations, 
rather than whole wave functions and technical issues of how they may be represented. 
And if applied carefully enough, they not only provide reliable self-consistent predictions 
but also link back to the full theory where they originate and thus provide fundamental 
insights. 

For many purposes, almost all the information contained in a wave function is irrele- 
vant, and a few state parameters of finite number suffice for all potentially observable con- 
sequences that can be imagined. This is what provides a much more economical derivation 
of physical results. One should note that effective equations are not merely an amendment 
of classical equations by quantum correction terms, although one can always obtain such 
equations in semi-classical regimes because the classical limit must be respected. However, 
effective equations apply more generally and constitute a systematic approximation scheme 
to analyze full quantum properties such as dynamical expectation values of a state. 

Best known among effective formulations are probably low-energy effective actions in 
particle physics [H [2]. But they are also available in canonical formulations [3j HI [5], 
where they have proven fruitful for quantum cosmology [6j [7]. They can be extended to 
constrained systems, where they provide effective constraints for the state parameters [§]. 
Thus, all ingredients are given which are necessary for an application to loop quantum 
gravity [9], [TUJ, [11] and a derivation of its effects. Due to the complexity of the problem, 
there is no complete derivation of effective equations for loop quantum gravity, but several 
characteristic quantum effects are known and can be analyzed. Taken together, all quantum 
corrections provide a complicated substitute for the classical equations, but they can be 
separated and studied individually. As we will see in this article, this provides crucial 
insights into what should happen in a consistent full theory of quantum gravity. We will 
explicitly construct anomaly-free constraints which incorporate quantum corrections of 
inverse metric components as they occur due to the discreteness of loop quantum gravity. 
A companion paper [12] will use these considerations of effective constraint algebras to 
provide quantum corrected cosmological perturbation equations in terms of gauge invariant 
variables. 

2 Effective equations and effective constraints 

We start by reviewing the scheme of effective equations for canonical formulations. For an 
unconstrained system, the quantum theory is given by a Hamiltonian operator H which is 
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self-adjoint on a given Hilbert space. It determines evolution of states if) by the Schrodinger 
equation 

ihip = H^ (1) 

and allows one to solve, e.g., for scattering amplitudes by evolving given initial states into 
possible final states. 

2.1 Effective equations of motion 

Alternatively, one can view the expectation value of the Hamiltonian, H® := (H), as a 
functional on the infinite dimensional space of states. It generates the same Schrodinger 
evolution by Hamiltonian equations of motion 

<<o H M. :{ <o>,^ } . (2) 

General expressions for the relevant Poisson brackets on the right hand side can be com- 
puted from commutators of basic operators used in the given quantum theory. 

A solution for (0)(t) to the Hamiltonian equation of motion has the same information as 
the expectation value (^\0\xjj) computed in a state satisfying the Schrodinger equation 
(or as the expectation value of a Heisenberg operator). However, in general H®, evaluated 
in a given state, depends not only on expectation values of basic operators but also on 
fluctuations and all other moments of the state. There is thus a complicated, infinite 
dimensional coupled system of equations involving not only the time dependence of (O) 
but also independent quantum variables such as (O n ) ^ (0) n (or expectation values of 
other operators) which appear in quantum theory. 

Effective equations are obtained when one can self-consistently determine regimes where 
these infinitely many equations can be decoupled to a finite system. This usually happens 
in semi-classical regimes where higher moments of a state are sub-dominant compared to 
low ones, but effective equations can be applied more generally. The 1-particle irreducible 
low-energy effective action, for instance, can be derived in an approximation consisting of a 
combination of an adiabatic expansion with one in h [3], [5]. Some rare solvable systems can 
be studied by exact effective equations without truncation, which in cosmology is realized 
for a flat isotropic model sourced by a free, massless scalar [6] . 

2.2 Effective constraints 

Gravity is governed by constraints rather than a true Hamiltonian. Just like the Hamil- 
tonian before, constraint operators Cj give rise to principal effective constraints (Cj), but 
with those an infinite tower of other constraints for quantum variables is generated [Sj. We 
are thus dealing with a system of infinitely many constraints on an infinite dimensional 
phase space even for a single classical canonical pair. The higher constraints can be ignored 
in our treatment of characteristic loop quantum gravity effects appearing in the principal 
constraints. As we will see, this is sufficient to arrive at a consistent constraint algebra 
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together with the cosmological perturbation equations it implies. Higher constraints would 
be required if we were interested in the constrained evolution of higher moments. 

The principal constraints are obtained from expectation values of quantum operators 
and thus contain several quantum effects. In general, they depend on quantum variables 
and include the quantum back- reaction of, e.g., fluctuations and correlations on expecta- 
tion values. But especially in loop quantum gravity they also contain characteristic effects 
which are a consequence of the fundamental quantum representation used. In loop quantum 
gravity, these translate consequences of the kinematical discreteness of the loop represen- 
tation into effective equations and thus show implications for dynamical states. The basic 
reason for properties of the loop representation is the use of SU(2)-valued holonomies of the 
Ashtekar-Barbero connection instead of linear functions of the connection |13j . In particu- 
lar, the basic variables become complex but still have to obey certain reality conditions to 
ensure that the correct number of classical degrees of freedom is quantized. This is usually 
implemented by requiring self-adjointness or unitarity of basic operators in the quantum 
representation, but one advantage of effective equations is that reality conditions can be 
represented representation independently at the level of expectation values and quantum 
variables [61 |8] . This remains true after solving the effective constraints, thus showing cru- 
cial properties of the physical Hilbert space even in cases where the physical inner product 
may be difficult to construct. 

A further advantage of effective equations, especially for the purpose of cosmological 
perturbation theory, is the issue of introducing a background geometry to define perturba- 
tive expansions in gravitational variables such as inhomogeneities. While the underlying 
theory of loop quantum gravity is background independent in a form which does not make 
it straightforward to introduce a perturbative background via states, effective constraints 
can easily be expanded by perturbing expectation values around a background of the de- 
sired classical form. The background then enters via a selection of a class of states to 
compute effective constraints [H]. (See also [15] for a conceptually similar proposal based 
on boundary states.) Using a Friedmann-Robertson- Walker background, for instance, al- 
lows one to derive cosmological perturbation equations directly — and effectively — from 
a background independent quantum theory of gravity. 

For the resulting set of equations to be meaningful, they must be consistent in that they 
derive from constraints which are anomaly-free also in the presence of quantum corrections. 
If this fails, it will be impossible to express the quantum corrected perturbation equations 
solely in terms of gauge-invariant variables as they are determined by the gauge flow 
of the corrected constraints. In particular, as we will demonstrate in this article, off- 
shell anomaly freedom is required. (The importance of the off-shell anomaly problem 
was also emphasized in [16] based on alternative fundamental considerations.) If one has 
an anomaly-free quantization with constraint operators such that the constraint algebra 
[Ci, C j] closes to a first class system, then also the algebra of principal quantum constraints 
{(Ci), (Cj)} will close because it derives from the commutator algebra. (If there are 
structure functions, higher constraints as mentioned in the beginning of this subsection will 
be involved.) Approximations to effective constraints then have to be done self-consistently 
in such a way that violations of closure of the algebra do not happen up to the order 
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considered in an expansion. 

In loop quantum gravity, however, no satisfactory form of all constraint operators is 
known which would satisfy the requirement of off-shell closure. (The arguments in [P71[T8] . 
for instance, specifically refer to partially on-shell statements, and also the reformulation 
of anomaly- freedom as a condition for the existence of observables in [191 EQ] is on-shell.) 
Without off-shell closure, on the other hand, physical applications based on the usual form 
of cosmological perturbation equations are impossible. (While applications may be possible 
based on a complete set of quantum observables in a form of reduced quantization, this 
route does not seem manageable.) The final advantage of effective formulations exploited in 
this article is then that one can ensure off-shell closure of the effective constrained system. 
Thus, one can include known quantum effects as they occur in a quantization where one 
has not yet taken care of anomaly freedom, obtain candidates for effective constraints 
with those corrections in a suitable parameterized form (reflecting either quantization 
ambiguities or incomplete knowledge of properties of a quantum operator), and compute 
their Poisson relations. In general, this algebra will exhibit anomalies, but in some cases 
one can adapt the correction functions used in the parametrization of quantum effects such 
that anomalies vanish. 

If there is no such adaptation, this specific quantum correction would be ruled out. But 
if one can successfully remove anomalies while keeping non-trivial quantum corrections, one 
will learn how specifically the quantum effect has to arise in quantum operators, and how 
completely quantization ambiguities can be fixed by the requirement of anomaly-freedom. 
The advantage of effective equations is then that one can do such an analysis order by order 
in various expansions, instead of having to face a complicated operator algebra in which 
all possible quantum effects are included at once. The result of completing such a program 
will not only be consistent sets of equations of motion which can be used for applications, 
but also valuable feedback on the underlying fundamental theory which in our case will be 
loop quantum gravity. Thus, we are providing a clear link between fundamental properties 
of quantum gravity and its phenomenology. 

2.3 Cosmological perturbation equations 

Linearization of Einstein's equations around Friedmann-Robert son- Walker (FRW) space- 
times provides cosmological perturbation equations for ten metric components. These 
metric perturbations are subject to coordinate (gauge) transformations parametrized by 
an infinitesimal 4- vector field £ M (/x = 0, . . . , 3), which, in presence of matter, generically 
give rise to six gauge invariant perturbations, i.e. combinations of metric and matter per- 
turbations which remain unchanged under linear changes of coordinates. In the linear 
regime, the former decouple into three independent modes: scalar, vector, and tensor, car- 
rying two degrees of freedom each. The evolution of vector and tensor modes taking into 
account corrections expected from loop quantum gravity was investigated in [21] and [22] 
respectively. In this paper, we focus on the scalar perturbations which, along with the 
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background FRW-metric, take the form 



ds 2 = a 2 (r]) (-(1 + 20)dr/ 2 + 2d a Bdrjdx a - ((1 - 2?P)5 ab + 2d a d b E)dx a dx b ) , (3) 

where the scale factor a is a function of the conformal time r\ and the spatial indices run 
from 1 to 3. The perturbations <p, ip, B and E are then combined into the two gauge 
invariant Bardeen potentials [23] 



# = ^-^[B-E} (4) 

$ = + (b - E) + (b - E 

whose evolution is governed by the linearized Einstein equations 

V 2 ^ - 3jr(jr$ + if) = -47rGa 2 5T 0{GI) (5) 
d a + 4f\ = -A7rGa 2 5T^ (6) 

-\d b d a (§ - *) = 4vrGa 2 5T a b(GI) . (7) 

Here a dot denotes derivative w.r.t conformal time, = - is the conformal Hubble 
parameter and 5T^ G1 ' are gauge invariant perturbations of the matter stress-energy tensor. 
These equations are commonly derived from the covariant field equations or by varying 
an action expanded to second order in the linearized fields. But Hamiltonian formulations 
exist for the same procedure, which is more suitable for a comparison with canonical 
quantum gravity (in particular in Ashtekar variables as used in [25]). 

To formulate the Hamiltonian setting, the action is used to determine Poisson brackets, 
and thus a decomposition into configuration fields and their momenta, as well as constraint 
functions. The constraints serve several purposes: (i) they restrict initial values of the fields 
to those allowed values which make the constraints vanish, (ii) they generate gauge trans- 
formations which in the case of general relativity agree with coordinate transformations, 
and (iii) they provide equations of motion for the fields in any coordinate time parameter. 
(The latter is itself subject to the coordinate changes by transformations generated by the 
constraints.) All this is necessary to reproduce a covariant system even though distinguish- 
ing momenta, which are related only to time but not space derivatives of fields, invariably 
removes manifest covariance from the Hamiltonian formalism. 

For this to be consistent, it is crucial that the constraints are preserved under the time 
evolution they generate. This is automatically guaranteed if they form a first class set, i.e. a 
set of functionals whose mutual Poisson brackets vanish when evaluated in fields satisfying 
the constraints. In other words, the gauge transformations and evolution generated by the 
constraints then define vector fields on field space which are tangent to the sub-manifold 
defined by the vanishing of constraints. Starting on the constraint surface, either changing 
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the gauge or following evolution will then keep us on the constraint surface. This is certainly 
realized classically, as a reflection of the general covariance of the underlying theory. 

However, if quantum aspects are implemented, one must ensure that this consistency 
requirement remains maintained: the quantization must be anomaly-free. Otherwise the 
equations may show the wrong type and number of degrees of freedom if formerly gauge 
quantities acquire gauge invariant meaning. Or, worse, anomalies may make the equations 
inconsistent to the degree that no non-trivial solution exists at all. Anomaly-freedom is 
thus a key requirement not only for the consistency of an underlying fundamental theory 
but also for the possibility of applications. Quantum corrections cannot appear in arbitrary 
forms, but only in restricted ways such that the constraints form a closed algebra under 
Poisson brackets. In particular, anomaly freedom will reduce some of the arbitrariness of 
the form of loop quantum gravity corrections. 

Moreover, as we will see explicitly, to provide quantum corrections to Eqs. ((SJ), © 
and ([7]) the algebra must close off-shell, i.e. it is not enough that the Poisson brackets 
of constraints vanish when the constraints are satisfied but even on parts of the phase 
space where constraints Cj do not vanish we must produce a closed algebra of a form 
{Ci, Cj} = fjj{A, E)Ck- (Here, A and E denote the canonical fields which may appear in 
the coefficients of the algebra; this means that in general we have structure functions rather 
than structure constants.) The effective algebra may differ from the classical one, and thus 
be quantum corrected as well as the constraints; but it must still close. The reason is 
that the whole set of coupled equations must be consistent, which presents a mixture of 
constraint equations ([5]) and ([6]) and evolution equations given by (J7J) together with the 
continuity or Klein-Gordon equation. To ensure that these equations are consistent, we 
must consider the constraints before they are solved. Consistency then requires an off-shell 
closure of the constraint algebra. Practically, the consequence is that only in this case 
we can express all the equations solely in terms of gauge-invariant variables as they are 
determined by the quantum corrected constraints. Once this is achieved, the equations 
are consistent and can be solved and analyzed. In the absence of off-shell closure, on 
the other hand, there would be left-over terms in the equations of motion which contain 
gauge-dependent quantities making such an evolution unphysical. 

In this context, it is important to realize that there is no shortcut to implementing the 
quantum corrections of fully perturbed field equations consistently. (Notwithstanding the 
fact that this has been attempted on numerous occasions such as [26j [27J [28] in the context 
of loop quantum cosmology, including by some of the present authorsjj) Consistency even 
for the purposes of phenomenological applications is intimately linked to the fundamental 
problem of anomaly-freedom once inhomogeneities enter the game. In homogeneous models 

particular, in Rcf. [28] a gauge fixing choice was made which, as follows from the present paper, 
resulted in anomalies. Also the recent proposal [29] of building inhomogeneous configurations from small 
numbers of patches subdividing space suffers from anomalies if more than two patches are considered. In 
this work, the authors also use a Born-Oppenheimer approximation — by assuming that the wave-function 
of the background metric and of the perturbations can be separated. This approximation is valid only if 
the wavelength of the perturbations A are much larger than the Planck length, i.e. A> lp. In the present 
work, such an approximation is not needed to derive the perturbation equations. 
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of gravity there is just one constraint, which clearly has a vanishing Poisson bracket with 
itself and thus forms an off-shell closed algebra. Thus, in homogeneous quantum cosmology 
there is no anomaly problem whatsoever. Here, quantum corrections can be implemented 
at will, only restricted by possible self-imposed conditions such as the desire to be as 
close to a candidate for a "full," non-symmetric theory as possible as it is expressed in 
loop quantum cosmology [30]. (Some of the structures, chiefly the kinematical quantum 
representation, of loop quantum cosmology can be linked to loop quantum gravity and 
are thus more restricted [HI (3TJ [321 [331 El]- But no such derivation exists yet for the 
constraints which are most important to see the precise role of quantum corrections on the 
dynamics.) 

It is then sometimes proposed to implement quantum corrections only in the background 
evolution, for instance by effects motivated from loop quantum gravity, and then use some 
inhomogeneous degree of freedom such as a matter field as a measure of perturbations 
around the background. If just the background is quantized, one knows corrections only in 
its evolution but would have to keep the structure of classical perturbation equations other- 
wise unchanged. This is rarely consistent, and the treatment is not gauge-invariant. Gauge 
invariant quantities in general relativity such as ([!]) combine several metric perturbations 
and possibly matter fields. Taking only a matter field, say, as the measure of perturbations 
means that one is fixing the gauge (by implicitly assuming non-gauge invariant metric 
perturbations to vanish) without even knowing what the gauge transformations are. The 
classical case of linear perturbations around Friedmann-Robertson- Walker spacetimes al- 
lows gauges where only the matter fields are inhomogeneous but not the metric like, for 
instance, the uniform density gauge. However, quantum corrections change the constraints 
and thus the gauge transformations they generate. The form and availability of certain 
gauges changes, and it is no longer possible to re-express the gauge-fixed results in terms 
of the gauge invariant quantities unless one considers the full gauge problem. This can 
only be done when initially all perturbations are allowed and the anomaly-issue is faced 
head-on. 

In the context of classical cosmological perturbations, the above arguments can be 
rephrased in the following manner. The effective corrections arising in loop quantum 
cosmology can formally be written as: 



where r^ iv contains all the corrections from loop quantum cosmology and T^ v corresponds 
to the stress-tensor of the classical matter field. Although the matter field might be an 
ideal fluid, the stress-tensor arising due the new physics cannot necessarily be treated 
as a perfect fluid. More importantly, the perturbation of the stress-tensor dr^ will in 
general have some anisotropic stress and the velocity perturbation 5v a will not vanish for 
a standard gauge choice. Hence, it is important to study the perturbations in a gauge 
invariant manner. 

If there is an anomaly-free version of quantum corrected constraints and the corre- 
sponding form of covariance, one could compute complete gauge-invariant quantities to 
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arbitrary orders in an expansion by inhomogeneities; see e.g. [351 EH [37] or, in deparam- 
eterized form after introducing dust as a clock matter system, [381 [39]. As discussed, this 
requires off-shell anomaly freedom of the constraints which is not easy to realize in closed 
form. The treatment by effective constraints then provides a key advantage: one can ver- 
ify anomaly-freedom order by order in the expansion by inhomogeneities (which may be 
combined with a semiclassicality expansion in H). This can be done with much more ease 
than a full anomaly analysis but still, as we will see explicitly, provides crucial feedback 
for the full theory. 

2.4 Correction functions 

Any quantization, such as loop quantum gravity, implies characteristic effects which change 
the classical behavior. Almost always, there are quantum back-reaction effects by state 
parameters such as fluctuations and correlations on expectation values. (The only excep- 
tions are free or solvable models such as the harmonic oscillator where moments of a state 
evolve independently of expectation values.) In addition, the specific quantum represen- 
tation may imply further characteristic effects, which in the case of loop quantum gravity 
are all related to the spatial discreteness of its kinematical representation. The classical 
set-up makes use of basic variables given by a densitized triad Ef, which provides the 
spatial metric via E^E\ = det(q c d)q ab , and the Ashtekar connection A l a = T l a + ^K^ with 



the spin connection and extrinsic curvature K\ = E^K a b/ y\det(Ej)\. This canonical 



pair of fields is then quantized in the form of fluxes, i.e. integrations of the triad over sur- 
faces, and holonomies or parallel transports of the connection. The resulting background 
independent representation has characteristic properties of spatial discreteness such as a 
discrete spectrum of flux operators (which contains zero). Such properties imply associated 
quantum corrections which appear whenever there are inverse powers of densitized triad 
components (some of which would classically diverge at singularities) or holonomies of a 
connection rather than just connection or curvature components. 

All these corrections typically occur at the same time and must be combined in a 
complete treatment. While one type of correction might be dominant in certain regimes, 
this would not be known a priori but had to be shown by a dedicated analysis. Nevertheless, 
due to the complexity of general quantum corrections, it is legitimate to separate the 
different corrections at first, analyze individual effects and then combine results. In spirit, 
this is similar to the calculation of corrections to an atomic spectrum, which can be done 
individually for relativistic corrections, spin-orbit interaction etc. and eventually combined 
in a complete spectrum. In this paper, we focus on loop quantum gravity corrections as 
they arise from inverse components of the triad. 

These corrections are already relevant for cosmology, where they have been analyzed 
in preliminary forms in homogeneous and inhomogeneous contexts pEOl HTl [26l W2\ H3~l HH 
SSI HH1 H7J HH] . (Note that a sub-dominance of these corrections compared to those due to 
holonomies has been claimed based on an analysis of isotropic models |49j. However, this 
is based largely on an inadvertent and artificial suppression in the models used [50]; see 
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also App. [A] In any case, the arguments put forward in the context of [49] do not apply 
to inhomogeneous situations.) The precise form of such corrections as they would result 
in a principal constraint from the expectation value of a constraint operator cannot, at 
present, be computed due to the complicated form of the volume spectrum which would 
be required (see e.g. [51]). But partially the behavior is known as it follows for instance for 
diagonal triads [52]. The typical behavior is that the classical function of triad components 
is multiplied by a correction function ot(Ef) which approaches the classical expectation one 
for large values of triad components. (More precisely, the function depends on fluxes, i.e. 
triad components integrated over elementary plaquettes of a discrete quantum state. This 
makes the functional behavior independent of the choice of coordinates.) 

At smaller scales, however, the function starts to deviate from one and implies quantum 
corrections. If the correction function is evaluated on an isotropic background [53], it has 
a peak at a certain characteristic scale a* of height larger than one, and then drops off at 
smaller scales to reach zero for vanishing triads. Notice that inhomogeneous contexts and 
states make it meaningful to speak about this behavior in terms of the scale factor a. In 
exactly isotropic models which are spatially flat, the absolute size of the scale factor has no 
meaning. However, the argument of correction functions is determined by a dimensionless 
ratio given by q := i\o? jl\ =: a 2 /a 2 where £q is the size in coordinates of an elementary 
plaquette whose flux appears as an argument. The product t§a has unambiguous meaning 
because it does not change under rescaling coordinates (which would change both a and 
£o individually). The peak of the correction functions occurs near q ~ 1, i.e. a ~ a*. The 
characteristic scale a* can be written as a* = £p / £ = (JV '/VoY^ip where J\f is the number 
of vertices of an underlying state contained in a region of coordinate volume Vq. The ratio 
N /Vq appearing in a* is thus the patch density of an underlying discrete state measured 
in a given coordinate system. For nearly homogeneous configurations, it does not depend 
on the region or on Vq, but on coordinates. (The physical vertex density which would be 
independent of coordinates is A/"/(a 3 Vo), but it would not be appropriate to determine a 
characteristic scale for a. Note that near a ~ a* there is one patch per Planck cube; upper 
bounds for the patch density can be obtained from phenomenological considerations, such 
as from big bang nucleosynthesis [54].) The value of a* depends on the normalization of 
the scale factor. But it also depends on the vertex density which can be large. Thus, the 
peak of correction functions for a denser state is realized on larger scales, which increases 
the corresponding quantum corrections. 

An additional implication of the appearance of the vertex density is that Af is typically 
history dependent [T4J E5] if the dynamical quantum evolution refines the state as the uni- 
verse grows (rather than just blowing up a fixed lattice). Thus, also the scale a* is history 
dependent which contributes to the regime dependence of this type of correction. For a 
given background, the history dependence can always be expressed as an a-dependence, 
which is sometimes seen as problematic because the scale factor is not coordinate indepen- 
dent. However, given that the origin of the refinement lies in the inhomogeneous setting, 
a proper reduction introduces the correct scaling dependence via additional parameters 
depending on the state; see also App. [A] 

For an implementation of perturbative inhomogeneities, regimes where relevant scales 
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fall below a* pose difficulties because the scale of inhomogeneity would be close to the 
discreteness scale. In this paper, we thus assume that scales of the densitized triad are 
above the characteristic scale a*, where correction functions deviate from one by terms 
perturbative in the Planck length: 



with positive coefficients c a and n a . Both coefficients can be derived from a specific quanti- 
zation but are subject to quantization ambiguities. The coefficient c a , in particular, is then 
related to a* (and to £ , providing the correct coordinate dependence in the presence of the 
scale factor). Thus, c a may itself depend on a if the vertex number Af in a fixed volume, 
and thus a*, changes with the universe expansion. We are assuming that the dominant 
a-dependence is via a power-law of the given form. 

Constraints for linearized perturbations will not only require the dependence of oi(Ef) 
on the triad when the latter is diagonal, but also the dependence on off-diagonal compo- 
nents. Classically one can always gauge the triad to be diagonal, but gauge transformations 
are quantum corrected and a consistency analysis of the equations must be done before 
a gauge is fixed. The off-diagonal dependence of a is not known in explicit form, and 
it is difficult to derive because unlike the diagonal case it requires non-Abelian features 
of the quantum theory [511 [56]. As we will see, the consistency analysis of constraints 
then relates the off-diagonal dependence to the diagonal dependence via the condition of 
anomaly-freedom. Moreover, other terms in the constraints, including matter Hamiltoni- 
ans, will also require characteristic quantum corrections which, in contrast to the primary 
correction functions, may not obviously be expected from explicit quantizations in homo- 
geneous models. Nevertheless, such additional corrections, called counterterms in what 
follows, are required for anomaly-freedom. In this way, they are fixed in terms of the 
primary correction functions depending on diagonal triads. All this not only provides 
consistent equations ready to be applied in cosmology, but also precise feedback on what 
terms a full anomaly-free quantum constraint must contain. As technical control on the 
full setting increases, these predictions will provide strong consistency checks of the whole 
framework. 

3 Canonical perturbation theory and primary correc- 
tion functions 

In this main part of the paper we develop the ingredients of a consistent perturbation 
theory in the presence of quantum corrections to the classical constraints. 

3.1 Constraints and primary corrections 

We first introduce the constraints and primary correction functions which are expected to 
arise in the effective constraints. Formally, the corrections are introduced as multiplicative 




(9) 
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factors of some terms in the constraints which depend on the phase space variables and 
approach unity in the classical limit. In this paper we restrict ourselves to correction 
functions resulting from the quantization of inverse-triad terms of the constraints. For 
the primary corrections, the functions are also assumed to depend only on the triad and 
to be local, i.e. independent of spatial derivatives of the triad. This reflects properties 
of these functions as they have been introduced in homogeneous models. The input can 
thus be used to formulate an initial expectation of the form of such functions. Moreover, 
in this section, we assume that the corrections can in principle be obtained from the full 
(non-perturbative) theory, and hence should depend only on the full triad Ef = Ef + SEf 
rather than on the background Ef and perturbations 5Ef as distinct arguments. Later 
on we shall analyze the consistency of such assumptions. Anomaly-freedom will generate 
additional counter-terms of further corrections, which can be re-interpreted as a connection 
dependence or non-locality of the primary corrections. Such a dependence is in any case 
expected for covariant corrections which can, e.g., be formulated as functionals of curvature 
invariants. Of course, we could put in such a dependence from the outset, but it would 
make the calculations much less tractable. 

General relativity in Ashtekar variables is subject to the Gauss, diffeomorphism and 
Hamiltonian constraints. The Gauss constraint is identically satisfied for scalar modes 
and does not need to be considered here. In the full quantum theory, the diffeomorphism 
constraint does not receive quantum corrections but the Hamiltonian constraint does [H] . 
The diffeomorphism constraint is thus taken as the classical one, D[N a ] = D graiV [N a ] + 
-Dmatter [N a ] with a gravitational part 

D gTav [N a ] := J d 3 xN a [{d a Ai - d b Ai)E) - A&Eft (10) 

and a matter part 

AnatterM = J d 3 xN a 7cd a¥ > (11) 

for a scalar field (p. 

We express the classical gravitational Hamiltonian as 

V[iV1 = i6^/ s d3xM (12) 

in terms of the Hamiltonian density 

n = ^== { - 2(1 + !- 2 )(A - r)*{A - r)f ) , (is) 

where the curvature of the Ashtekar connection is given by 

F k ab = 2d [a A\ + Zii k KA{ , 

7 is the Barbero-Immirzi parameter and the spin connection r„ is considered as a functional 
of the densitized triad (written explicitly in Eq. (128]) ). The presence of an inverse of the 
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triad determinant, whose quantization does not have a direct inverse because it has a 
discrete spectrum containing zero, suggests the presence of a primary correction function 
of inverse-triad type, 

H^ V [N] = -|— / d 3 xNa(E?)H = H^ y [aN] =: tf grav [A>] . (14) 

For the same reason, primary quantum corrections v{Ef) and cr(Ef) are introduced into 
the matter part of the Hamiltonian constraint as 



H 

matter 



[N] = J d 3 xN {vH* + aU v + H v ) , (15) 



where 



^Wf' ^^^> < 16 > 

are the (classical) Hamiltonian densities corresponding to the kinetic, gradient and poten- 
tial terms respectively. Note that the potential term is not expected to acquire primary 
quantum corrections because there is no inverse of the triad in this term. Thus, the cor- 
rection does not simply amount to a rescaling of the lapse function even if v and a would 
equal a. 

In App. |B]we discuss the Poisson brackets between unperturbed H P [N] = H p av [N] + 
^matter W) and D[N a ] = D zrav [N a } + letter [iV a ], as well as between H^N,} and H P [N 2 ]. 



3.2 Perturbations 

Consider first an action which depends only on one scalar field (p. Generalizations to 
arbitrary tensor fields will be considered in the following subsection and in App. After 
the Legendre transform the action takes the form 

S[(p] = Jd 4 x(7T<f-H{<f,7r)), (17) 

where TC is the Hamiltonian density and tt is the field momentum. Given a space-time 
slicing by constant-time surfaces, we split the fields (p and tt into their homogeneous parts 

<p := y Jd 3 x <p, tt := y Jd 3 x vr (18) 

and the inhomogeneous remainder 

S(p :— if — if, 5tt := 7i — tt. (19) 

Here, Vq is the volume of a spatial slice if it is closed, or can be thought of as a very large 
(but finite) infrared cutoff volume otherwise. The coordinate size Vq will only appear in 
basic variables and their symplectic structure, but not in final equations of motion. 
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We also require the inhomogeneities 5ip and 5n to be small: 



5<p 



<1, 



Sir 

TT 



< 1 



(20) 



for the slicing we use, so that this can be considered as perturbations around homogeneous 
solutions. For a generic Hamiltonian, these conditions may at some point be violated during 
the evolution. In fact, only a narrow class of Hamiltonians admits such a splitting, for which 
the inhomogeneities remain smaller than the mean fields at all times. At the moment, as is 
common in cosmology of the early universe, we merely assume that (1201) holds for the regime 
under consideration. Hence, from now on we shall refer to 5ip and 5n as perturbations and 
will also speak of the first, second, etc. perturbative order, denoted by superscripts (1), 
(2), . . . in what follows. Specifically, we will be interested in the perturbations up to the 
second order in the Hamiltonian, which implies a linear perturbation theory in terms of 
equations of motion. 

From the very definition of 5tp and 5tt it follows that any first-order quantity averages 
to zero. In particular, 



Xi := / d^xX^tp = 0, xi ■= / d 3 x\ 2 5n = 



(21) 



where Ai and A2 are 'smearing' constants!! Therefore the first term in the action ( TT7T) splits 
into two parts: 



J d 4 x7T0 = j d 4 a;(7r + dir)(<p + dip) = V J dtinp + J d 4 x5rr5ip, 



(22) 



yielding the basic Poisson brackets 



{ < P,k} = —, {6<p{x),6ir{y)} = 5 3 {x - y) 



(23) 



and, for phase space functions, {, } := {, }^ + {, }stp,Sni where 
{F, G}^ 
{F,G} s <p,6k I (V s ./- 



1 fdFdG _ dFdG 
Vq \d(p dfc dit dip 
5G 



, ( 5F 5G 5F 8G \ 
" \6{6<p) 6{6ir) ~ 5{5n) 5{5<p) ) 



(24) 



As discussed in App. [Cj these brackets are not fully general and in some cases care may 
be required, but they are sufficient for calculations done here. 



2 Obviously, the constant factors Ai and A2 are not necessary in the constraints. They have been 
introduced to make a more explicit connection with the case of tensor fields, where such factors would be 
necessary to contract spatial and internal indices 
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3.3 Perturbed constraints 

So far we have used the Ashtekar connection A\ as one of the canonical variables as 
required for holonomies. From now on we will explicitly use the spin connection and 
extrinsic curvature, 

A i =V + iK l 

where 7 is the Barbero-Immirzi parameter. Also the canonical pair K % a = K l a + 5K l a , E? = 
Ef + SEf can be split into the homogeneous parts 

K = k8i, E«=p8l 

corresponding to the flat FRW-background, and the inhomogeneous perturbations which, 
for the scalar mode, are described by a pair of scalar functions each: 

6K* a = <5>i + d a d i K 2 , 6E? = Sfex + did a e 2 . (25) 

Note that in the perturbed context, the independent phase space variables are (k,p) and 
(5K l a ,5E^), and the non-trivial Poisson brackets between them are given byB 

{kp} = ^, {5K{x),5E b ] {y)}=^G5)5 b a 5\x-y) . (26) 

The Hamiltonian density (fT5|) , expressed in terms of the extrinsic curvature, becomes 

n = 4 k |pt= [2d c r d + 4»(iTr3 - ktk)] + n 7 , (27) 

y|det£/| 

where the last term is proportional to the Barbero-Immirzi parameter, 

n = 7e f 2E ! E * ( dK i + 4 n r^jQ) = 7e f 2E l E t D K i = 2 g D gi 

7 1 v^det^ V d m " c ^ 4 v^det^ d v^et£j 

with the Gauss constraint (? J , which thus vanishes. Indeed, the Gauss constraint implies 
that the extrinsic curvature can be written as K l d = K^E^ / yj det E\ where K dc = K c d- 
Consequently, 

rpc rpd rpb 

7 3 det£ 

Thus the classical theory in (K l a , Ej) is explicitly insensitive to the Barbero-Immirzi pa- 
rameter, as it should. The 7-dependence, however, will appear in the correction functions 
resulting from a quantization procedure (after which no unitary transformation exists to 
change 7 without leaving a trace on observable quantities). 



3 See also App[C]for details. 
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The remaining part of the Hamiltonian density can be expanded straightforwardly in a 
perturbation series, although the spin-connection requires some care. Its full expression is 

K = ~\e- k E) (d a E$ - d b E* + EtE&El - , (28) 

where E\ with a lower spatial index designates a co-triad of density weight minus one, 
whose perturbed expression reads 

K = -$1 - -^SE%5 ca 5 kl . 
p p z 

The first order part of (1281) . 

Sr a = 1 (ep b a - e*5l + e^8 ac + d b 5E° , (29) 

is simplified significantly for a scalar perturbation of the form ff25l) . The diagonal part of 
5Ej (in the term dbSEj, which is SjdbEi,) contributes to the linearized spin connection 

ST? ias) = 1 (0 - ejftei - e?te + 3e?ftei) = iejftei . 

For the off-diagonal perturbation, on the other hand, the expression d b 5Ej = d b d c djS2 is 
symmetric in the indices b, c,j, implying that only one term in (1291) remains: 

Jr n(ofr-diag) = 1 _ o + o + e »^Ae 2 ] = ^e%As 2 . 
2p 2p 

Combining the last two expressions, we obtain the linearized spin connection 

5r a = 8rJ> d iag) + srf ff - diag) = — Sd, ( e i + A ^) • ( 30 ) 

2p 

Note that this expression is diffeomorphism invariant to linear order. Remarkably, the 
(gradient of the) term in the parenthesis can be expressed as the divergence of the unsplit 
triad perturbation 

dj (e 1 + Ae 2 ) = d a 5E«, 

which can be easily checked by inspection. Thus, for scalar mode the linearized spin 
connection can be expressed as 

5r a = ^d b 5E]. (31) 

The second order part of the gravitational Hamiltonian constraint also contains a term 
quadratic in r". However, as such a term is necessarily multiplied with a background 
quantity, the term becomes proportional to the trace of the spin connection, 5*1^. For 
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the scalar perturbation, the latter can be shown to vanish up to at least third order using 
similar symmetry arguments. 

In the spin connection part of the Hamiltonian, the first order term is contributed solely 
by the 'derivative term' 



2e 



rpc rpd 

jk j k d c sr d 



y/\ detE\ 



(i) 



Vp 



&d a 8E° 



whereas the second order part comes from both the derivative and the quadratic terms: 



c rpd 



2e 



jk E j E k 
; vTdetE 



S c 5Y d 



(2) 



V 



3/2 



and 



jk j k j r rm r rn 



y/\ detE\ 



(2) 



1 



2 p3/2 



S ij d a 5E^d b 5E b i . 



Combining the last two terms, we obtain (up to a total divergence) 

(2) 



e* S5L (2d c 5T d + e\JT™5T n d ) 



yf\ detE\ 



1 



2j]i/d 6»d a 6E«d b 5E». 



(32) 



Expanding also the extrinsic curvature term, we thus arrive at the expression for the 
gravitational Hamiltonian density H = + TiP^ + with 



-6k VP , 



I 2 

y/p" 



Vp5Ki c 8K k d 5 c k 8f - VPiSK^) 2 - -=6E?6Ki 



2k 

w 

-ic rj\2 



(33) 



Likewise, the perturbed diffeomorphism constraint including up to quadratic order in 
perturbations is 



-Dgrav [N 



8ttG 



d 3 x5N c [pd c (5 d k 5K k d ) - p(d k 5K k c ) - U k c {d d 5E d k )] . (34) 



We now consider the contribution from the scalar matter sector. The classical Hamil- 
tonian is given by 

(35) 



^matter [N] = J d*xN (H„ + H V + H v ) , 
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where the kinetic, gradient and potential terms are defined in (|T6|) . Again, we have a 
perturbation expansion with 



«f'=^-^^. < = °. «?>=^(V,(^ + ^^) (37) 



and 



* ~ 2 p 3 / 2 p 3 / 2 2p 2 p 3 / 2 \ 8p 2 Ap 2 



= Vv w (^)V+p 3/ v,(^)^ c J 



The perturbed diffeomorphism constraint for the scalar matter field is 

^matter [iV a ] = J d 3 x6 N c 7rd c 6<p . (39) 

In the following two subsections we explicitly compute the Poisson brackets between the 
perturbed classical constraints and show that their algebra is closed. At the same time, 
we will be including primary correction functions and see how the algebra changes. 

3.4 Poisson bracket between Hamiltonian and diffeomorphism 
constraints 

We begin by considering the gravitational sector of the classical constraint algebra. We will 
see later that for computational purposes it is convenient to split the classical perturbed 
gravitational Hamiltonian as 

# g rav[iV] = -±- / d 3 xNH = H svav [6N] + H giav [N] . (40) 
lbnG ./v 

Here if grav [£iV] includes only the perturbed component of the lapse function whereas 
Hgra-vlN] involves only the background lapse. Explicit expressions for each part of the 
perturbed Hamiltonian constraint are 

H S ^ V [N] = -J— [ d 3 xN [W<°> + H {2) ] , H giav [5N] = -A- / d 3 x5NH m , (41) 

1D7TG J 107rG J 
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where perturbed Hamiltonian densities are given in equations (1331) . We consider now the 
Poisson bracket between the gravitational Hamiltonian H grav [N] in (j4"Tl) and the gravita- 
tional diffeomorphism constraint D grav [N a ] in ( J34l) : 

{H grav [N] , D grav [N a ] } = -H grav [5N a d a 5N] . (42) 



This Poisson bracket between classical perturbed constraints (|42|) is very similar to its 
counterpart between the full classical constraints [TO], also computed in App. [B] This 
demonstrates the consistency of perturbed constraint expressions and elementary Poisson 
brackets between background and perturbed basic variables. 

As in the gravitational sector, the classical perturbed Hamiltonian for the scalar matter 
field including up to quadratic terms in perturbations can be expressed as H majtter [N] = 

^matter [SN] + H mattev [N] where 

^matter [5N] 1= j d?x5N [?#> + W«] . (43) 

Perturbed Hamiltonian densities for scalar matter are given in equations fl36|) . (1371) and 
( |38l) . The Poisson bracket between the matter Hamiltonian constraint and the total dif- 
feomorphism constraint can be computed as 

{# matter [N] , Dgrav [N a ] + Matter [N*] } = "Matter [6N a d a 6N] . (44) 

Combining gravitational sector ( 1421) and matter sector (14"4"1) contributions, we can evaluate 
the Poisson bracket between the total Hamiltonian and diffeomorphism constraints as 

{H[N],D[N a }} = {H gTav [N], D grav [N a ]} + { # matter [N], D grav [N a ] + D matter [N a ] } 

= -H gTav [5N a d a 5N]-H mattei [5N a d a 5N] = - H[5N a d a 5N] . (45) 

Clearly, perturbed expressions of total constraints along with elementary Poisson brackets 
between background and perturbed basic variables satisfy the same Poisson brackets as 
the full expressions. 

We now analyze the situation for primary corrected constraints. As in the classical 
situation, we split the primary quantum corrected gravitational Hamiltonian constraint as 

H^ V [N] = Y ^Jd 3 xNa(Et)n = H^ av [5N] + H^ v [N}- (46) 

The part of Hamiltonian constraint containing only the perturbed lapse function H giav [SN] 
and the part of Hamiltonian constraint containing only background lapse H grav [N] are 
defined as 

H*[N] := ^/dy^ + a^^ + a^ + a^) 
s 167rG J 

H^JSN] := ^ J d'xSN + a^H®) , (47) 
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where a = a^°\ The Poisson bracket between the primary quantum corrected Hamiltonian 
constraint (|46l) and the diffeomorphism constraint (1341) can be computed as 



{H*[N\, D pKf [N a ]} = -HF{6N a d a 6N]+A i 



HD 

grav 



(4S 



where 



grav 



16ttG 



d s x{d c SN j )p 



5NH {0) da(1) (5 a 5 c - 5 c 5 a ) + N Ih {0) ( ' -— 



1 da 
_l 

3 9p p 



+ 



da® 



d(5E: 



5*5?: 



+ U {1) — — (5 a 5 c 

d(6E?y 3t 



(49) 



would appear as an anomaly if primary corrected constraints were used as quantum con- 
straints: the Poisson bracket ( 1481) between the quantum corrected Hamiltonian constraint 
and diffeomorphism constraint has additional terms which cannot be expressed completely 
in terms of the gravitational constraints for any lapse function or shift vector. Thus, these 
terms in the constraint algebra are potentially anomalous. 

Next we explore this issue for the quantum corrected scalar matter sector. Similarly to 
the classical Hamiltonian constraint, the quantum corrected matter Hamiltonian can be 
split as 



matter 



[N] 



J d 3 xN [u(E9)7U + <r{E*)H v + H v ] =: H^ tt J5N] + H, 



p 

matter 



[N] .(50) 



The two parts H^ atteT [SN] and H^ atter [N] of the matter Hamiltonian are defined as 



H, 



matter 



[N] 



d s xN 



{mf + nf) + [vM-h!® + i/Mww + + m% ] + n 



H: 



matter 



[5N] 



d 3 x5N [vWHW + mp + nW] 



(51) 



where v = u^°\ Here z/W, and denote zeroth, first and second order terms 
in perturbations of the quantum correction function v. The Poisson bracket between the 
quantum corrected scalar matter Hamiltonian floDl) and the total diffeomorphism constraint 

matter 

[iV a ] can be computed as 

{Fitter [N] , ^grav [N*] + D matter [N a ] } = -Fitter [5N a d a SN] + A™ teI . (52) 

As in the gravitational sector, there are additional terms present also in the matter sector 
Poisson bracket which are 



A 



HD 



- / d?x(d c 5N j )p 



bNH^^— {8*51 - 8«5f) + N ( 5 C - 



d(5E? 



P 



+ 



i dv m dv® 



3W + wW) m ~ m J + ^ WW) m ~ ^ 



(53) 
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Matter sector anomaly terms are similar to anomaly terms in the gravitational sector, 
but there are important differences. In particular, matter anomaly terms involve only the 
kinetic sector of the matter Hamiltonian density TC n . In contrast, gravitational anomaly 
terms contain the total gravitational Hamiltonian density TC. Thus, one cannot even hope 
to combine all anomaly terms to form the total Hamiltonian constraint for specific cor- 
rection functions. Moreover, cancellation would require a lapse depending on correction 
functions. The requirement of an anomaly-free constraint algebra then demands that 
gravitational sector and matter sector anomaly terms must vanish separately. Combining 
contributions from the gravitational sector ( f4"8l) and the matter sector ( 1521) . we can write 
the Poisson bracket between the quantum corrected total Hamiltonian constraint and the 
total diffeomorphism constraint as 

{H P [N], D[N a ]} = —H p [5N a d a 5N] + <f v + A^ ttev . (54) 



3.5 Poisson bracket between two Hamiltonian constraints 



We now consider the Poisson bracket between two Hamiltonian constraints smeared with 
different lapse functions. It can be split into three components as follows 



{H grav [JVi] , # grav [N 2 ] } + {// ma ttcr [#i] , # matt er [N 2 ] } 
+ [{^matter [JV X ] , H gIav [N 2 }} - (N, <- N 2 )} . (55) 



Using the perturbed expression of the classical gravitational Hamiltonian (f4T!) we compute 
the Poisson bracket between gravitational Hamiltonian constraints as 



{tfgravt^^W^]} 



{H grav [SNt] , H grav [N] } + { H grav [N] , # grav [SN 2 ] } 
{H^ISN, - 5N 2 ], H grav [N]} = D grav jd a {5N 2 - 5N X 



(56) 



where we have used the property that {H gTSuV [5Ni\, H grav [5N2]} = 0. Similarly, using the 
perturbed expression of the classical scalar matter Hamiltonian ( )43i) . we compute the pure 
matter sector contribution as 



{H mat ter[Nl}, H mattev [N 2 ]} = D 



— d a (6N 2 -6Ni 
V 



(57) 



It is easy to show that the net contribution from the Poisson bracket between gravitational 
Hamiltonian and matter Hamiltonian parts in the constraint vanishes. In particular, 



{fUttofJVxLiWiVa]} - (^i ^ iV 2 ) = . 



(5* 



Combining equations fl56|) . fl57|) and fl58l) we evaluate the Poisson bracket between total 
Hamiltonian constraints as 



{H[N 1 ],H[N 2 ]} = D 



N 
P 



d a {SN 2 - 6Nx) 



(59) 
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Thus, the perturbed expression of the classical Hamiltonian constraint indeed satisfies the 
same Poisson bracket with itself as its unperturbed expression. 

Now using the perturbed expression of the primary quantum corrected Hamiltonian 
(146p we compute the Poisson bracket between quantum corrected gravitational Hamiltonian 
constraints as 



{H^m, H^ av [N 2 }} = {H^N.-SN^H^N}} = D 



grav 



a 2 — d a (6N 2 -6Ni 
V 



+A 



HH 

grav 

(60) 



with 



A 



HH 

grav 



8nG 



d 3 xAT(5AT 1 



5N 2 



^rw. > \ + + 3^) 



dp 3 d(5Ef> 1 



(da da^ 



da® 



■ (61) 



Similarly, using equation (I50|) we compute the Poisson bracket between primary quantum 
corrected scalar matter Hamiltonians as 



{-^matter [^l] ^matter [^2]} = D 



ua—d a (SN 2 -5N 1 ) 
V 



(62) 



Eq. f[62|) is analogous to its classical counterpart except that the new shift vector for 
the resulting diffeomorphism constraint now contains quantum correction functions v and 
a. Net contributions from the Poisson bracket between quantum corrected gravitational 
Hamiltonian and matter Hamiltonian constraints are 



No 



A 



HH 



(63) 



with 



A 



HH 

matter 



tfxNiSN! - 5N 2 



TV 2 



7T 



P 



3/2 



(2kVp) t- 



dv 



2p 3 /2 



<9p s(<^; 



7T 



2^3/2 



(2/CVP) 



d(6E? 



6? 



+ 



7T 



A' 



4<9z> 



# + 



3 dp c d(5E s 



d{5Ef 
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(64) 



2^3/2 ^/^ 

Using equations (160]) . fl62l) and fl63|) we can combine contributions from the gravitational 
and matter sectors to express the Poisson bracket between primary quantum corrected 
total Hamiltonians: 



{H^N^H^N,}} = D 

-\- -Dmattcr 



N 



a 2 —d a (5N 2 - SNi) 
V 



{ud -a 2 )—d a {5N 2 -5N 1 ) 
V 



+ AHH , JHH /gKN 
T /V„ av T -^matter ■ l UJ / 



grav 



22 



3.6 Conditions for an anomaly-free constraint algebra 

In contrast to the classical situation, we have seen that primary quantum corrected con- 
straints fail to form a first class constraint algebra for arbitrary correction functions. To 
interpret this properly, we recall that quantum correction functions that we have used as a 
guideline in the Hamiltonian constraint are not completely known. In particular, one can 
compute only zeroth order terms using homogeneous models [531 E3 EE] (as well as some 
partially gauge-fixed inhomogeneous cases using lattice states [52]). Linear and quadratic 
terms in perturbations of quantum correction functions can in principle be computed using 
the machinery of the full theory [51]. However such computations are not yet available. 
In this section we will analyze whether there are conditions on quantum correction func- 
tions that we must impose based solely on the requirement of an anomaly-free constraint 
algebra. 

First, we note that the explicit appearance of the matter diffeomorphism constraint in 
equation (1651) . drops out if the quantum correction functions satisfy 

a 2 = vo . (66) 

This requirement may be seen as a consistency relation between gravitational and matter 
correction functions, which has important physical implications: For instance, it ensures 
that gravitational waves and massless scalar fields propagate with the same group velocity 
given by the physical speed of light |22j . 

Furthermore, in order for the constraint algebra to be closed, Ag T ^ v should vanish ir- 
respective of the choice of lapse function. In other words, anomaly terms involving back- 
ground lapse N and perturbed lapse 5N must vanish independently. This requirement 
leads to 

On the other hand, from equation ( IBTi) the conditions 

Ida <9a (1) <9a (1) <9a (2) 

3^ a = WW) and ~df = WEff (68) 

are required to make A^ v vanish. Clearly, the requirement of an anomaly-free constraint 
algebra imposes restrictions on the first and second order terms of the quantum correction 
functions. However, it is evident that equations (I6T|) and (|68|) are over-complete for the 
unknown functions a' 1 ' and a^ 2 \ Importantly, these conditions are incompatible with each 
other for non-trivial primary corrections and admit only trivial solutions of a = constant, 
a^ = and = which is just the classical situation without quantum corrections. 
The situation for the scalar matter sector is very similar. Using equation (1531) it is easy to 
see that A^ v = requires 

ldi?6E° «9i/ 2 ) 
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while 

sW i ~ WW) ~W ~ WW) 1 ' m 

solves ^matter = 0- As In the gravitational sector, anomaly-free requirements on matter 
quantum correction functions admit only trivial solutions of v — constant, v^> = and 
z/( 2 ) = 0. There are no additional requirements on the quantum correction function o as 
only its background component appears in the perturbed Hamiltonian. However, equation 
( |66l) then requires that also a must be constant. 

At this stage, we could only conclude that inverse triad corrections leave no trace what- 
soever in effective constraints of an anomaly-free quantization. This would be extremely 
puzzling given the crucial role played by the corresponding operators for well-defined funda- 
mental constraints of loop quantum gravity. Fortunately, what we have shown is, in fact, a 
weaker statement since we assumed the primary correction function to depend only locally 
on the triad in algebraic form. What we have shown is that this pure triad dependence is 
insufficient, and we will now relax this assumption by introducing additional corrections 
which we call counterterms. These terms are not directly motivated by simple expressions 
computed from a constraint operator, but they will be fixed in terms of primary correction 
functions by anomaly freedom. In the conclusions, we will comment on the expectations 
for the presence of such terms based on a loop quantization. 



4 Anomaly-free quantum constraints 

In the previous section, we have shown that the presence of only primary quantum correc- 
tion functions does not lead to an anomaly-free perturbed constraint algebra. It is however 
possible that the chosen form of quantum correction functions, as they have been used in 
all studies so far, does not capture all possible quantum effects. Naturally, one is then led 
to ask whether there are "counter terms" to the chosen form of the correction functions 
that should be included in quantum corrected expressions of the Hamiltonian constraint 
to make the constraint algebra anomaly-free. We show here that such expectations are 
indeed realized. In particular, it is possible to arrive at a quantum corrected constraint 
algebra which is anomaly-free by including specific counter terms to the primary quantum 
corrected Hamiltonian constraint. 



4.1 Gravitational sector 

For a non-trivial primary quantum correction function it is not possible for both A^ v and 
^grat to vanish simultaneously. However, it turns out that one can perform partial anomaly 
cancellation in the constraint algebra even for non-trivial quantum correction functions by 
relaxing some of the conditions imposed so far to result in manageable computations. 
We approach this by adding counter terms, i.e. further potential quantum corrections, 
ensuring first that the quantum corrected Hamiltonian constraint is covariant under spatial 
diffeomorphisms, i.e. A^ v = 0. At this point, it is important that the diffeomorphism 
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constraint should not receive quantum corrections of the type studied here. Counter terms 
thus appear only in the Hamiltonian constraint. We have seen that condition (1671) on the 
quantum correction function, in particular cr 1 ) = 0, precisely ensures this requirement and 
we can simplify the primary gravitational constraint (45) to 

K^m ■■= d 3 ^Ar [awW] . (71) 

Similarly, we simplify the expression of (59), which we then refer to as 

p p p ^ 

As we have already seen, canceling this anomaly based solely on primary corrections 
could be achieved only for the trivial case of constant a. To cancel the anomaly terms (1721) 
without having to require constant a, we need to generate additional terms which take a 
form similar to primary anomaly terms. In fact, additional corrections not considered so 
far can easily arise in an effective Hamiltonian constraint. Adding additional terms into 
the Hamiltonian constraint however can potentially generate new anomaly terms in the 
Poisson bracket {H^. av [N}, _D grav [iV a ]}. To avoid this we should ensure that the counter 
terms being added to the Hamiltonian constraint commute with the diffeomorphism con- 
straint. Moreover, counter terms should not affect the background dynamics as there are no 
anomalies in the constraint algebra when one turns off inhomogeneity. Thus, counter terms 
should be constructed using only those terms which contain inhomogeneous perturbations. 

We notice that while we describe the emergence of counter terms in a constructive 
manner, all this reflects requirements which fundamental anomaly-freedom would pose. To 
that end, we consider a minimal approach for constructing counter terms. In particular, 
we require that counter terms should generate only those three kinds of terms which are 
already present in the anomaly expression (1721) . With this requirement, the allowed form 
of counter terms that can be included in the quantum corrected Hamiltonian is given by 
HCJN)=HCJ5N]+HC &V [N) where 



-Af(p)k^p(8]8Ki) - g(p)^(6l6E< 



and 



(73) 



(74) 



Here we have introduced three dimensionless scalar functions /, g and h which depend on 
the quantum correction functions and are to be determined through anomaly cancellation 
conditions. Only background components of these functions are relevant as the counter 
terms are already quadratic in perturbations. 
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The new terms can be interpreted as resulting from a dependence of the primary a on 
extrinsic curvature components and spatial derivatives of the triad as a general functional. 
Thus, the introduction of counter terms relaxes some of the conditions imposed earlier on 
a. The Poisson bracket between counter terms and the diffeomorphism constraint can be 
computed as 

1 



{H^[N\, D[N a ]} 



ItiG 



d 3 xa(d c 6N c )6N [-(2f + g)(k 2 ^p 



(75) 



Thus, the requirement that counter terms commute with the diffeomorphism constraint 
leads to the simple condition 

2f + g = (76) 

on the coefficient functions. Out of three unknown coefficient functions only two remain to 
be determined. Using equations (1751) . (173|) we compute contributions from counter terms 
to the Poisson bracket between Hamiltonian constraints 
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(77) 



Since the perturbed classical constraint algebra is closed without counter terms, the 
quantum counter terms must vanish in the classical limit, i.e. when all primary correction 
function are unity. Counter terms can then depend on primary correction functions only 
through their derivatives. Given the expanded form (Q of primary quantum correction 
functions used here, terms such as (da/ dp) 2 can be neglected compared to the terms 
da j dp. For the same reason the contributions from the Poisson bracket between counter 
terms {-ff^ av [iVi], H^. av [N 2 ] } can be neglected compared to the contributions considered in 



Combining all contributions from the original anomaly 
(177)1 one can express the total gravitational anomaly ^4 grav : : 



21) and from counter terms 
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Anomaly cancellation will require coefficients of (d c d^5Ej), (SjSK^.) and (5 J c 5Ej) to vanish. 
This in turn implies that the coefficient functions /, g and h should be such that they 
satisfy three equations Q\ = 0, Q2 = and Q3 = 0. On the other hand, / and g also 
need to satisfy equation ( 1761) . Thus, the set of equations for /, g and h may appear to be 
over-complete. However, it is remarkable to note that equation (1761) along with equation 
flHUj) solves the equation (181]) identically. In particular using equation (1751) . it is easy to see 
that 

g 3 = -Q2 ■ (82) 

Thus for a given quantum correction function a, there are unambiguous solutions for /, 
g and h such that the constraint algebra is anomaly-free. In particular, for a background 
quantum correction function a given in 

, = 9 = 2 P da h = 1 2n Q + 1 p da 
2 An a + 3 a dp ' 4n Q + 3 a dp 

solve equation (1761) and ensure vanishing of equations (1791) . (180]) and (18I1) . 



4.2 Cosmological constant 

For obtaining anomaly freedom in the gravitational sector by adding appropriate counter 
terms, it was crucial that coefficients Q2 and Qz are related through equation ( 1821) . Thus, 
one should consider the robustness of this relation under the inclusion of additional classical 
terms. Including a non-zero cosmological constant to the gravitational sector provides a 
definite test to see whether such a relation can still be satisfied. Counter terms in the 
gravitational sector would now generate additional contributions due to the presence of 
the cosmological constant term. We now show that a non-zero cosmological constant does 
not spoil the non-trivial consistency condition (182]) . 

Contributions to the gravitational Hamiltonian constraint from a non-zero cosmological 
constant A are 

H a[ N \ = -TTi f d 3 xNy/\detE\A =: H A [5N] + H A [N] . (84) 

07TG J 

As with the matter potential term, no primary inverse triad corrections are expected. 
The perturbed expressions of H\[5N] and H\[N], including up to quadratic terms in 
perturbations are given by 

Ha[N] = ^ J d 3 xN [nf } + Hf] , H A [5N] = ^ I d 3 x5NH A 1} . (85) 
The explicit expressions of perturbed densities 7~C£ , l~l A and Ti A are 

Hf ) = 2Ap 3 / 2 , ^ = 2Ap 3 / 2 ^-S (86) 

(2) , , /2 /(^^) 2 (5 k c 5 J d 5E^Ei)\ t 

n m = 2Af/ 2 _ ico (g7) 



27 



Contributions to the anomaly expression arising from the Poisson bracket between counter 
terms and cosmological constant terms are 



' ' d 3 xiV(<5Aq - 5N 2 )a{Ap) 



8nG 



-(micbf-k 



2p u yj 



Combining contributions from counterterms and the original anomaly in the presence of a 
non-zero cosmological constant, one can evaluate the total gravitational anomaly ^4g rav := 



-^•grav -^grav A&. as 
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where new coefficients in the anomaly expression are 
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19) 



(90) 



(The A-dependence cancels upon using the background Friedmann equation.) Using equa- 
tion (1761) . which remains unchanged, we again note that the new coefficients satisfy the 
same non-trivial consistency relation 



(91) 



Thus, also in the presence of a non-zero cosmological constant there are unambiguous 
solutions for /, g and h such that the constraint algebra is anomaly-free. This demonstrates 
that anomaly freedom of the quantum corrected constraint algebra including appropriate 
counter terms is a robust feature. 



4.3 Scalar matter 

For cosmological applications, we must ensure the existence of consistent equations in 
the presence of matter. Similarly to the gravitational sector we ensure first that the 
quantum corrected scalar matter Hamiltonian is covariant under spatial diffeomorphism, 
i.e. ^matter = 0- This requirement implies that we should impose conditions f l69l) on 
the quantum correction function v and simplify the expression of the primary quantum 
corrected matter Hamiltonian as 

(pn® + nf) + (z/W«jw + m m + + m <$ + nf) 

HL^N] := Jd 3 x6N[iMV+HU] . (92) 



HLuJN) := /d 3 xiV 
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This also simplifies the matter anomaly term A^ tter which we then refer to as 
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(93) 



For computational convenience we consider the construction of counter terms for the kinetic 
and potential sectors separately. 



4.3.1 Kinetic sector 

To cancel anomalies in the kinetic sector of scalar matter, we start with a general form of 
possible counter terms as H° V [N] := H%[5N] + H$[5N] + H%[N] + H$[N) where 



H%[5N] = / d 3 x5N 
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g 5 (p)^5 ab d a 5cpd b 5cp 



(95) 



The general guidelines followed for the gravitational sector led us to introduce seven di- 
mensionless unknown functions fi, / 2 , gi, g 2 , 93, g^ and g 5 which should be related to 
primary quantum correction functions as it will be determined through anomaly cancella- 
tion conditions. 

The Poisson bracket between counter terms (j9"4l) . (|95|) and the total diffeomorphism 
constraint is 



{H^[N],D[N a ]}= I d*x 



(d c 6N')N g 4 — 



( z/TT 2 6E?\ f U7T 2 

^w-^) + i8JNC) \ SN{2fl - h) w^ 



+N( 9l - 92)—^ - N(2g 2 -g 3 + g A ) 



utt 2 (5i5E" 



P 



2 p3/2 2 p 



.(96) 



Requiring that counter terms commute with the diffeomorphism constraint leads to the 
conditions 

2/i = h , 9i = 92 , 2^2 = fl's , gi = . (97) 

Given that g 4 is required to vanish, we will drop the corresponding term from the set 
of counter terms in our further evaluation. We began with seven unknown functions in 
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the counter terms for the kinetic sector. Requiring that counter terms commute with the 
diffeomorphism constraint imposes four conditions. This in turn allows only three more 
conditions to be imposed on these functions from anomaly cancellation in the remaining 
Poisson bracket of Hamiltonian constraints. 

There are some subtleties in finding anomaly cancellation conditions for the matter sec- 
tor. Inclusion of counter terms to the gravitational sector has generated additional contri- 
butions both in the gravitational as well as the matter sector. Thus, matter sector counter 
terms need to cancel the original anomaly expression (193]) but also contributions from 
gravitational counter terms. Contributions from gravitational counter terms to anomaly 
expressions of the matter kinetic sector are 
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Similarly, contributions from counter terms of the matter kinetic sector to the Poisson 
bracket between total Hamiltonians can be computed as 



A c 



{H%[N x l H'[N 2 ] + H^[N 2 ]} - (Ni 



- v 



'grav 



d 3 xiV(5JVi - 5N 2 



-ffV 2 ^(/! + g 5 ) 
P 



N 2 ) 
vtiStx ak 



+ 2p3/ 2lj c) Vp [ H) 2p 3 / 2 2p VP 



pi I a ^/p y Up 

-(-2p^ + 4f 2 -3g 3 



(99) 



We then combine the original anomaly (|93|) with contributions (|98|) from gravitational 
counter terms and contributions from matter kinetic sector counter terms (|99|) to express 
the total anomaly A^ := ^matter + ^S-avTrV + i n the kinetic sector of scalar matter as 
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In the presence of a non-zero scalar matter potential, the imposition of background and 
perturbed Hamiltonian constraints does not determine matter kinetic terms in terms of 
gravitational terms in the Hamiltonian constraint. In such situations anomalies in the ki- 
netic sector of scalar matter must vanish independently of the gravitational sector anomaly. 
From, Eq. (HOOD it is evident that anomaly freedom in the kinetic sector requires four con- 
ditions to be satisfied, i.e. B\ = 0, B 2 = 0, £> 3 = and £> 4 = 0. However as we mentioned, 
after imposing equations fl97j) we have only three undetermined functions in kinetic sector 
counter terms. Thus it may appear once again that there is over-determination of counter 
terms. However, similarly to the gravitational sector, coefficients of the anomaly expression 
(1 1001) satisfy a non-trivial consistency relation. In particular, using relations ( 1761) . ( 1971) . one 
notes that 

B 3 = -2B 1 - B 2 . (105) 

Thus, cancellation of kinetic sector anomalies requires only three equations to be satisfied 
by counter term coefficients. In other words, counter terms for kinetic sector are unam- 
biguously determined by anomaly cancellation conditions. 



4.3.2 Potential sector 

We recall that the potential sector did not involve any primary quantum correction function 
and did not contribute to the matter sector anomaly. However, including counter terms 
in the gravitational and matter kinetic sectors leads to new anomaly terms involving the 
scalar matter potential. Such new anomaly contributions from gravitational counter terms 
to the Poisson bracket between Hamiltonians can be computed as 
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Similar anomaly contributions from counter terms of the matter kinetic sector are 
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(107) 



Thus, counter terms of the matter kinetic sector generate a new anomaly term involving 
Vipipi'fi)- Gravitational counter terms on the other hand, do not lead to such a term. 
For non- vanishing f\, not all terms can cancel by combining equations f ll06p and (11071) . 
Even though we did not consider primary quantum corrections in the potential sector, for 
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anomaly freedom we need to allow counter terms even here. As in the kinetic sector, we 
begin with a general expression of possible counter terms in the potential sector H^[N] := 
H°[5N\ +H$[N] where 



H?[6N] = / d 3 x5N 
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We have introduced six new unknown functions fa, fa, g%, g-j, g$ and g$ in the counter 
terms of the potential sector. To ensure invariance of the potential sector counter terms 
under diffeomorphism constraint we compute the Poisson bracket between counter terms 
and the diffeomorphism constraint: 
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It is then easy to see that counter terms commute with the diffeomorphism constraint only 
if the coefficients satisfy 



fa = , 5-7 = 0, 5-8 = 0, 59 = 



(111) 



Thus diffeomorphism invariance of the counter terms allows just two independent functions 
in the potential sector. The non-vanishing counter terms in the potential sector reduce to 
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Contributions from counter terms of the potential sector to the Poisson bracket between 
total Hamiltonians can be computed as 
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Combining (11061) . fjlOTj) and (11 1 3 P we form the total anomaly term A v := A^^+A^+A^, 
in the scalar matter potential sector: 
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As in the gravitational and matter kinetic sectors, coefficients of the potential sector 
anomaly satisfy a non-trivial consistency relation 



-Vo 



;ii6) 



using fl97j) . We recall that counter terms f II 1 3 j) of the potential sector have only two 
unknown functions ^3 and g§ which would be determined by choosing, say, T>\ — and T> 2 = 
0. Then, we would automatically have V3 = 0, but there are still non-vanishing terms in 
the anomaly fll 14j) of the potential sector. In particular, last two terms in fll 14j) are similar 
to the anomaly terms due to cosmological constant (jSSJ) and can be absorbed into anomaly 
terms of the gravitational and matter kinetic sectors by subtracting the total Hamiltonian 
constraint with suitable lapse function from it. Vanishing of total anomaly terms from 
potential sector then requires that should also vanish i.e. T>± = 00. However, requiring 
X>4 to vanish imposes additional restriction on counter terms which in turn requires primary 
correction functions a, v and o to satisfy another consistency requirement (see Eq. fll34j) ) 
in presence of a non-zero scalar potential, apart from the relation (1661) . After imposing 
V>\ = 0, T>2 = 0, T>3 = and P4 = 0, we can express remaining terms in equation (11141) as 
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where part of the potential sector anomaly to be included in the gravitational sector 
anomaly can be expressed as 
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4 The anomaly term involving D4 can be absorbed into anomaly terms of the gravitational and matter 
kinetic sectors by suitably subtracting total Hamiltonian constraint from it. However, such subtraction 
does not lead to a consistent set of solutions for anomaly cancellation conditions. 
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and another part that needs to be included in the anomaly expression of the matter kinetic 
sector is 

2 ak ,r;r (8i5E(, 
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We have seen earlier that the presence of a non-zero cosmological constant modifies the 
anomaly cancellation conditions as reflected in equation fl90|) . Similarly, the presence of a 
non-trivial scalar matter potential leads to changes in the anomaly cancellation conditions 
for both the gravitational sector as well as the kinetic sector of matter. In particular, we 
can combine equations f!78|) and (11181) to express the total gravitational anomaly *4. gra v : = 
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Using equation (J75]l . we again note that the new coefficients satisfy the same non-trivial 
consistency relation 

(122) 



Thus also in the presence of a non-trivial scalar matter potential, there are unambiguous 
solutions for /, g and h such that the gravitational sector of constraint algebra is anomaly- 
free. Similarly, for the matter kinetic sector we can combine the original anomaly (1931) . 
contributions (I98j) from gravitational counter terms, contributions (I99l) from kinetic sec- 
tor counter terms and contributions (11191) from the potential sector to express the total 
anomaly in the kinetic sector of scalar matter as 
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2^3/2 2p 

with new coefficients 

Bf = B 1 , B% = B 2 + f 

in the matter kinetic sector anomaly expression. Also here, the new coefficients satisfy a 
consistency relation 

Bl = -213% - B% (125) 

using (J75|) and flDTl) . This relation is analogous to equation (I105p and likewise it ensures 
that anomaly cancellation conditions lead to unambiguous expressions for counter terms 
in the kinetic sector of scalar matter when there is a non-trivial potential. 
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4.4 Quantum corrected total Hamiltonian constraint 



We now combine all primary correction functions and the counter terms to form the quan- 
tum corrected total Hamiltonian constraint H®[N] := H P [N] + H C [N] for the system 
consisting of a scalar matter field with arbitrary potential: 

H Q[N] := -J— / d 3 xiV [«H«°» + a^H Q ^ + aH Q ^] , 
s 1o7tG J 

H^ V [SN] := ^ J d 3 xSN [aH Q M] , (126) 

where the background density is unchanged except for the explicit factor of the primary 
correction, i.e. TC^^ = . However the perturbed densities now involve counter terms: 
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H Q W = ^Ki5K k d 5^ - VP(5K15 C 3 ) 2 - ~^>^K (127) 
/<2 ;5E<5E d k 5% + -^{5EW c f - (1 + h)£L(dJ5EZ)(d d 8Et) . 



It should be noted here that the terms in gravitational Hamiltonian which involve coun- 
terterms, contain either trace or divergence of perturbed basic variables. Thus, inclusion 
of counterterms does not affect the earlier results for vector and tensor modes [21, 22J. The 
complete quantum corrected matter Hamiltonian constraint is given by 

H^ ttcI [N] = jf d*xN [(^(0) + ^Q(o)) + (> W Q(o) + m m + - n m + h q {2) j 

H^uJSN] = J d 3 x5N [m^ + , (128) 

where background densities are again unchanged, i.e. = Ti^ and = . As 

in the gravitational sector, perturbed matter densities include counter terms and are given 
by 

«?" - d + A^-d + A)^^ (129) 
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?£3(2) = pi/2 



v 



To summarize the conditions on non-vanishing coefficients of the counter terms, we 
note that there are three such functions in the gravitational sector (11271) . six in the kinetic 
sector and two in the potential sector of scalar matter (11291) . Thus for the system under 
consideration we have a total of eleven correction functions contained in all counter terms. 
Invariance of counter terms under diffeomorphisms, (1761) and (|97|) ■ led to four conditions 

g = -2f , f 2 = 2h , g 2 = 9l , g 3 = 2g 2 (130) 

among the non- vanishing coefficients. These equations trivially lead to the solutions for g, 
f 2 , g 2 and g$, leaving seven functions to be determined. 

Cancellation of anomaly terms from the Poisson bracket between Hamiltonian con- 
straints led to three conditions (I12ip from the gravitational sector 

Qf = , Ql = , Qt = , (131) 

among which only two are independent due to (I122p . These two independent equations 
explicitly solve / and h in terms of the primary correction function a. Thus, there are only 
five remaining functions that need to be determined. Anomaly cancellation from matter 
kinetic sector (11241) leads to four conditions 

Bf = , B% = , B% = , B\ = . (132) 

Given the relation (I125p . there are only three independent equations which leads to explicit 
solutions for fi, gi and g§ in terms of primary correction functions a and v. The remaining 
two free functions f 3 and g^ are constrained by requiring cancellation of anomalies from 
the potential sector which gives four conditions 

V l = , V 2 = , £> 3 = , £> 4 = . (133) 

Using equation (11151) . one notes that D\ = and V 2 = already determines both / 3 and 
#6 in terms of other counter terms coefficients which are already fixed. 

While D 3 = is not an independent equation, X> 4 = imposes a non-trivial restriction 
on counter terms. Since all of them have been determined at this stage, consistency requires 
the primary correction functions to satisfy 

a p p ( a p\ up p I v ' p\ 2p I v p^ 



a 3 V a J v 9\v J 9 \ v J K J 

Note that this relation ties the matter correction function to the gravitational correction 
function, but it is independent of the matter fields. Finally, from (l66l) we have the relation 
a 2 = va to be satisfied by the primary correction functions. 

To summarize, the requirement of anomaly freedom in the constraint algebra tightly 
controls the allowed forms of primary correction functions. For primary corrections of the 
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form (Q, for instance, one can easily see that solutions exist provided certain relations 
between the coefficients c a , c v and powers n a , n u in the two primary correction functions a 
and v are satisfied. Thus, quantization ambiguities are non-trivially reduced, which allows 
stringent consistency tests by direct calculations from a full representation of the under- 
lying operators. These restrictions indirectly help to eliminate some of the quantization 
ambiguities encountered in quantizing inverse triad operators. 

5 Conclusions 

In this paper we have analyzed quantum corrected constraints at the perturbative effective 
level. The key issue has been whether they form a closed Poisson algebra, which would 
ensure consistency of equations of motion they generate. We have found that a behavior 
of correction functions a, v and a (scalars of zero density weight) as in homogeneous 
models, which would suggest that they (i) depend only on the triad Ef (but not on the 
extrinsic curvature Kl), (ii) depend on the triad algebraically (i.e. do not contain spatial 
derivatives of Ef), and (iii) in the perturbed context, depend on the background triad Ef 
and its perturbation 5Ef only in the combination Ef + 5Ef = Ef, would allow a closed 
algebra only if the corrections are trivial. 

At the same time, from the constructive point of view, it is not surprising that the 
three conditions cannot be met together. Indeed, the only scalar quantity that can be 
constructed from the triad alone is its determinant — a density weight one object — leaving 
no possibility of cancelling the density weight. One could relax any of these conditions and 
see whether that would allow non-trivial corrections. For instance, by allowing correction 
functions to depend on spatial curvature, which would require spatial derivatives of the 
triad, we could alleviate the problem of zero density weight since this would make available 
the quantity 

Ef Q ( d b E\ \ 
VI detE | a \ ^\AetE\) ' 

On the other hand, if we relax the first condition, quantities of the form EfK % a /^\ det E\ 
would be allowed. 

We were therefore led to conclude that expectations from homogeneous models did not 
capture all possible quantum effects, and turned to investigating what quantum corrections 
of inhomogeneous constraints would be allowed in principle, and which ones should be ruled 
out. In this process, we have generated several counterterms in addition to the primary 
corrections suggested by homogeneous models. The resulting counterterms admit non- 
trivial quantum corrections, and their presence and form can be related to fundamental 
aspects of loop quantum gravity. For instance, we have seen that quantum corrections 
must be connection dependent even when they come from inverse triad corrections. This 
can be interpreted as meaning that the computation of effective constraints, based on 
expectation values of constraint operators, must be done in coherent states such that a 
holonomy dependence of inverse triad expressions results. Correction functions must also 
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depend on spatial derivatives of the triad, which can be seen as leading terms in a derivative 
expansion of non-local expressions involving fluxes, i.e. 2-dimensionally integrated triads. 
There are also quite unexpected effects, such as counterterms in the matter sector involving 
derivatives of the potential. Not all of them are simply realized as a consequence of the 
expansion of V((p + Sip) by inhomogeneities. This suggests that the matter potential must 
be quantized in a non-local way to ensure anomaly freedom. This form of non-locality 
is currently not realized in quantizations of scalar matter in loop quantum gravity. It 
suggests concrete ways to change full constructions so as to provide an (off-shell) anomaly- 
free formulation. 

From the perturbed second order constraints one can directly derive Hamiltonian equa- 
tions of motion for the perturbed variables as well as gauge transformations on them. 
Both equations of motion and gauge transformations will be corrected by quantum gravity 
terms, which has to be combined for equations of motion of gauge invariant variables of the 
form (ISTl-fTj). Imposing the conditions found for an anomaly-free constraint algebra must, 
on general grounds, result in a consistent set of equations. This has already been verified 
for vector and tensor modes (see [21] and [22] respectively). In a companion paper [12] we 
explicitly derive gauge transformations and construct gauge invariant variables taking into 
account quantum corrections, which we will then use to derive gauge invariant equations 
of motion describing cosmological perturbations. 

We have provided one consistent set of equations by a process which demonstrates 
that possibilities of non-trivial quantum corrections are rather tight. In fact, existing pro- 
posals for primary correction functions are non-trivially restricted. Yet, different versions 
may be available, which could in principle be compared with full derivations of effective 
Hamiltonians to fix remaining ambiguities. But there may also be quantization ambiguities 
which cannot be removed based solely on consistency considerations; they would have to 
be restricted phenomenologically instead. It is thus important also for a fundamental un- 
derstanding to evaluate cosmological implications of the quantum corrected perturbation 
equations. 

In addition to other choices regarding one type of quantum corrections, which in this 
paper is inverse triad corrections, there are different general types of corrections. In loop 
quantum gravity, we have two additional classes: corrections of higher powers of the con- 
nection or extrinsic curvature due to the use of holonomies, and genuine quantum back- 
reaction effects which include the influence of the whole wave function on its expectation 
values. (It is the latter which underlies constructions such as the low-energy effective ac- 
tion used in particle physics.) These corrections turn out to be more difficult to compute 
in consistent form, which is still in progress. Our consistent equations are thus not to be 
considered as complete effective equations, and including the corresponding terms of one 
type may add to the effects of another type or decrease them (in a way which is regime 
dependent). But it is unlikely that complete cancellations happen because corrections of 
the different types take so different forms. Moreover, a complete cancellation would mean 
that the characteristic fundamental representation of loop quantum gravity would leave no 
trace on the physics of the theory. 

While quantitative results are expected to depend on the specific form of corrections and 
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the interplay of different types, the occurrence of qualitative effects signalling deviations 
from classical relativity is more robust. This differs from other results of loop quantum 
cosmology, such as bounces in homogeneous models where a sharp zero-result for the time 
derivative of a scale factor is required. Such sharp conditions can easily be destroyed when 
additional quantum corrections are included; see e.g. [7j. Compared to that, the complete 
elimination of qualitative effects of one type of correction by including another type is 
highly unlikely. 

The consistent constraint algebra shows that non-trivial quantum corrections which 
reflect the underlying discreteness of spatial geometry are possible. In this sense, general 
covariance is preserved. However, we have shown that the classical constraint algebra, 
while consistently deformed, is not represented exactly but receives quantum corrections 
from the corrected constraints. One can see this directly from (!60|) . for instance, which 
carries a factor of a 2 in the smearing function of the diffeomorphism constraint. This 
is required by consistency since the classical algebra cannot be realized with non-trivial 
quantum corrections. Thus, an effective action of loop quantum gravity cannot be sim- 
ply of higher-curvature type. (Non-local features, for instance, would then be essential.) 
Nevertheless, we expect that some of the corrections can be formulated by effective higher- 
curvature actions which applies even to the inverse triad corrections used here. Some of the 
counterterms, which depend on extrinsic curvature components as well as spatial deriva- 
tives of the triad, can in fact be interpreted as bringing the corrected constraints in a form 
amenable to being formulated as a higher-curvature action. In this context, we emphasize 
that the absence of new degrees of freedom in this Hamiltonian framework is not in conflict 
with the higher-derivative nature of higher-curvature effective actions as also discussed in 
[5j: a perturbative interpretation of higher-derivative actions, which is the only appropriate 
way in quantum gravity, does not give rise to more solutions than expected classically [55] • 

Acknowledgements 

This work was supported in part by NSF grants PHY0653127 and PHY0456913. MK was 
supported by NSF grant PHY0114375. SS wishes to thank Roy Maartens and Kevin Van- 
dersloot for discussions. He is being supported by the Marie Curie Incoming International 
Grant IIF-2006-039205. 

Appendix 

A Comparison with isotropic models 

An important measure for the size of quantum corrections is the characteristic scale a* 
which signals the onset of non-perturbative effects. As a critical value for the scale factor, 
it does not have absolute meaning because it can be rescaled by a choice of coordinates. 
It is ratios such as a/a* which have physical meaning related to the patch density of a 
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quantum gravity state. For a denser state features of correction functions based on inverse 
triad components are realized on larger scales, which increases the corresponding quantum 
corrections. 

These effects are also important if one tries to include the behavior in homogeneous 
models, even though an exactly homogeneous model provides only limited means of re- 
ferring to spatial discreteness of an underlying state and its refinement. For this reason, 
care is needed if one tries to address possibilities of refinement schemes and the size of 
quantum corrections in purely homogeneous settings, as is often done due to the simplicity 
of homogeneous models. Quantum corrections in a fully inhomogeneous situation must be 
expected to be larger than in a naive isotropic quantization which ignores the factor M 
of the patch density and implicitly assumes A/" 1 / 3 ~ 1 as in [60]. This is the reason why 
some minisuperspace considerations artificially suppress those corrections. Corrections 
from holonomies, on the other hand, increase with decreasing vertex density such that 
they would appear to be more pronounced. It is possible to mimic the enhancement of 
inverse triad corrections even in exact homogeneous models by computing their correction 
function for operators based on higher representations of SU(2) instead of the fundamental 
one [57J. The corresponding spin label is then related to the vertex density. 

In addition to the size of corrections, there is also the issue of the correct scaling 
behavior of correction terms. To have independence of the coordinate size Vq of the region 
whose patches are counted, we must have Af oc Vq. However, this provides coordinate 
independent quantum corrections only if we multiply by another function which can absorb 
the coordinate dependence of Vq. The simplest possibility in isotropic models is to use 
M oc a 3 Vo for which corrections depend neither on the size of the volume nor on coordinates. 
This behavior is indeed well-motivated based on lattice refinements (where the physical 
vertex density is constant) and was introduced in [51] based on scaling arguments. If N is 
allowed to depend directly only on the scale factor (and not on a/ a, say), and there is no 
other parameter which rescales under changing coordinates other than Vq, this is indeed 
the only consistent choice. In this sense, the behavior proposed in [61] is unique in isotropic 
models up to a single constant which determines the absolute number of patches. 

This uniqueness is, however, contingent on conditions which are too strong for reliably 
modelling what happens in inhomogeneous situations. If the model is no longer exactly 
homogeneous, the refinement of underlying states is history dependent in ways which do 
not simply amount to a dependence on a. One can always express the refinement as a 
dependence on scaling- independent observables such as a /a which provide an equally good 
measure for the history of different phases of the universe. For a given background solution, 
one can then express this as a dependence on a alone, given that all observables depend on 
a. In general, however, this provides more complicated functions than just A/"(a) oc a 3 Vo 
for the patch density. Some phases can, for instance, be described by a power-law form 
J\f(a) = Aoa^Vo where Ao arises in a complicated process by expressing the refinement via 
a function only of a. In particular, because the original refinement is history dependent 
only via observable quantities, this constant will automatically be equipped with a scaling 
dependence such that Nq(l x Vq is coordinate independent even if x ^ 3. The emergence 
of such a parameter can only be seen in the proper inhomogeneous context, invalidating 
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considerations based solely on homogeneous models. 

B Poisson brackets between unperturbed constraints. 

It is instructive to compute Poisson brackets between primary quantum corrected con- 
straints without expanding by inhomogeneities. We first consider the gravitational Hamil- 
tonian constraint (JHj). As H^ av commutes with matter diffeomorphism -D ma tter> ^ * s suffi- 
cient to compute {H^ TlLV [N], D grav [N a ]}. Classically, this Poisson bracket is 

{# grav [iV], D giav [N a }} = -H giav [N a d a N] . (135) 

The (gravitational) diffeomorphism constraint acts as a Lie derivative on a (gravita- 
tional) phase space function 

{F(A,E),D grav [N a }} = C Na F . (136) 
The Poisson bracket between Hg TSuV [N] and L>gr av [A^ a ] then is 

{H^ av [N],D giav {N a ]} = {[ d 3 xNaH,D giav [N a ]} = [ d 3 xNjC N «(aH) 

Jt Jt 

= [ d 3 xN (N a d a (aH) + {d a N a )aH) = [ d 3 x{-N a d a N)aH 
Jt Jt 

ee -H^ v [N a d a N}. (137) 

Where we have used the fact that the quantity &H has density weight one to expand the 
Lie derivative and integrated by parts in the next line. By the same token, we would not 
obtain the correct algebra if a would not be of density weight zero. Clearly, any functional 
of the (gravitational) variables defined by integration must be an integral of a density- 
weight-one function. Thus the only restriction on the correction function, obtained so far, 
is that it must be of zero density weight. In fact, only this condition makes the spatial 
integrals well-defined. 

The matter Hamiltonian constraint has non-zero Poisson brackets with both gravita- 
tional and matter parts of the diffeomorphism constraint. However, the total diffeomor- 
phism constraint acts as a Lie derivative on any function of all phase space variables 

{F(A, E, <p, 7r), D[N a }} = C N aF . (138) 

Hence its Poisson bracket with the (gravitational) diffeomorphism constraint should boil 
down to an expression analogous to ( 11371) . Again, the only condition on the correction 
functions is that they have zero density weight. In that case not only are the quantum 
corrected constraints first class, but also form an algebra identical to the classical one. 
In what follows we will make extensive use of 
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Lemma 1 Consider a functional 

F[N] = I d^x N(x)f((p,n) (139) 



of two canonically conjugate scala^ fields cp and it. If f does not depend on spatial deriva- 
tives of the fields, the Poisson bracket 

{F[iVx],F[iV 2 ]} fe7r) = (140) 

vanishes 

Proof: Since the integrand does not contain spatial derivatives, we have the functional 
derivative SF[N]/Scp = Ndf /dcp and 

□ 



On the contrary, if spatial derivatives of the fields are present in the integrand, the 
relevant functional derivative involves derivatives of the smearing function which implies 
a non-vanishing final expression for the Poisson bracket after anti-symmetrization over N\ 
and N 2 . 

Using this result, let us analyze the expression 

{H^Nd, H P [N 2 }} = {H^JN,} + H^JN^H^JN,] + H^ tt jN 2 }} 

= {H^MH^JN,}} + {/Ctter^J, Ct.^} ( 141 ) 

+ ({if£ av m , H^ ttcI [N 2 ] } - (N, <- N 2 )) 
term by term. The gravitational constraints yield 

= ^ grav [iV 1 9 iV 2 -iV 2 9 a iV 1 ]+^ gravgrav , (142) 

where we have used the fact that if N± t 2 were independent of phase space variables then 
we would simply have the classical constraint algebra but with new lapse functions Ni and 

5 Generalization to the case of tensorial as well as several conjugate fields is straightforward. 
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N 2 . However, since Ni j2 do depend on the densitized triad there is an extra (potentially 
anomalous) term in the Poisson bracket which is the second term ^4g rav grav, proportional 
to the derivatives of the correction function. The non-trivial contributions to the anomaly 

-W = - / = d'x {n^H± ( I ^3*^) - W ~ N.) 

- \ NA ( aiv S) - {Ni - *>} (143) 

come from the gradient terms of the Hamiltonians. Note that, for convenience, we switched 
the order of terms in the first line of (j!43p . In the second line, the only term in the 
parenthesis that survives the anti-symmetrization is the one proportional to the gradient 
of the lapse function d c N2. Thus the anomaly simply boils down to 

Aravgrav = H? av [M a ] with M a :=2e ijk -^=J=-— ; (N 1 d c N 2 -N 2 d c N 1 ). (144) 

V det E o^k 

It is easy to see that the symmetricity condition 

da da 

E °3Ei = tHE} (145) 

is sufficient to make the anomaly (j!44j) vanish due the contraction with e^. We should 
point out that (11451) is definitely satisfied for any triad-dependent scalar function, which 
has all internal indices contracted. 

The cross Poisson bracket {if^ av [iVi], H^ autteT [N 2 \} — (Ni <-> N 2 ) can be computed sim- 
ilarly. In the absence of curvature couplings, such that the matter Hamiltonian contains 
neither connection nor spatial derivatives of the triad, this Poisson bracket is given by 



No 



3 | £#grav[iVl] da 

(1 ■' ' SAi dE? N2Hv 



H^[M v } + H^[M a } 




where 



H?[M U \ = f d 3 xM„Hir, H?[M a ] = [ d 3 xM a Hv (147) 

JT, JT, 

with the effective lapse functions 



(i4s; 
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similar to (11441) . These vanish if the correction functions satisfy 



dv 



pa dv 

k dEf 



E* 



da 
3 9Ef 



F a da 



(149) 



Finally, the Poisson bracket between two matter Hamiltonians involves only functional 
derivatives with respect to the matter variables cp and it. By virtue of lemma ([I]), the 
non-trivial contribution comes from 



dFxNu/H*, / d 3 xN 2 aHv > - (Ni 



No 



Since the correction functions do not depend on the matter variables, they act as constant 
factors, i.e. 

{#ltter[iVl], tfmattert^]} = Anatter [v(j{N 1 d a N 2 - ^9°^)], (150) 



Combining (11421) . (11461) and (11501) and assuming (I145j) and fll49j) we obtain 

{^[iVx], # p [iV 2 ]} = D grav [a 2 (iV 1( 9 a iV 2 - A^A^)] + Anatter[^(iVicPiV 2 - jV^JVOKlSl) 

It is easy to see that the constraint algebra closes, if a 2 = va in addition to the requirement 
that a, v and a are all scalars of vanishing density weight. In that case the right hand side 
of (I15ip reduces to the total diffeomorphism constraint 

{H P {N X \ H P [N 2 }} = D[a 2 (N 1 d a N 2 - N 2 d a N l )] = DfiVi^iVa - N 2 &*Ni)]. (152) 

So far in this appendix, we have worked non-perturbatively which gives only a few 
conditions on quantum correction functions. The anomaly freedom conditions (|67|) and 
(!68|) obtained in the main part of this paper, where the condition of vanishing density 
weight turns out to be quite non-trivial, appear much more restrictive compared with the 
relatively mild-looking requirement derived in the context of the unperturbed system in 
this appendix. It is therefore pertinent to comment on this apparent discrepancy. 

Note that the conditions on the three correction functions imply the same functional 
form of a, v and a. Thus we shall restrict our consideration to only one of them. In section 
13. 1| we had made the following assumptions concerning the primary correction function a: 

(i) a depends only on the triad Ef (but not on the extrinsic curvature K\ or the connection), 

(ii) a depends only algebraically on the triad E" 1 (but not its spatial derivatives) 

(iii) in the perturbed context, a depends on the background triad E°- and its perturbation 
8Ef only in the combination Ef + 5Ef = Ef (i.e. a is expected to originate from a full 
unperturbed expression). 

One can check by inspection that assumption (iii) implies that (!68l) is automatically 
satisfied. Indeed, using the Taylor expansion 



a(Ef) = a(Ef) + 



da(Ef) 



dEf 



5E? + 



1 d 2 a(Ef ) 



2 dEfdE\ 

1 J 



5Ef5Ef+--- = a+a^ + 



(153) 
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it is easy to see that the terms on the right hand side are not entirely independent. Clearly 
the relations between a, a^ l \ and a^ are exactly written in Eq. (1681) . 

However, of greater concern is the other condition, Eq. (|67j) . In particular, the re- 
quirement a^ 1 ' = (along with (1681) ) rules out all possible non-trivial solutions. In order 
to understand its origin let us revisit the seemingly trivial restriction on the correction 
function to be of zero density weight. We start by formulating the following 

Lemma 2 A scalar ce(E?) of density weight zero satisfying the three assumptions above 
must be a constant function. 

Proof: Consider a scalar a(E£) of density weight w satisfying the aforementioned assump- 
tions. Its Lie derivative along an arbitrary shift vector N a is given by 

£-a = N b d b a + wad b N b . 

On the other hand, 

C^a = ^C*E? = ^ (N%E? - E\d h N« + E«d b N b ) . 

These equations are valid for any iV a . In the context of cosmological perturbation theory, 
N a = 0, hence N a = 5N a . In the perturbative expansion of the right hand side of the 
equations there is no contribution from the background part. Equating the corresponding 
linear order terms, we obtain 

( da \ (0) 

wad b 6N b = ( — (E?d b 6N b - E b d b 6N a ) . 

Using Ef = p5f the derivative (da/dE?)^ = (da/ f dEf)\ga can be rewritten as ^6fda/p, 
which yields 

^l^ db6N b = W(X d b 5N b . 
3 op 

The divergence of a generic shift vector does not vanish, and therefore the derivative of the 
background correction function is da /dp = ^pwa. Requiring w = results in da /dp = 
and consequently from (|68|) . a^ = 0, a^ = and so on. This concludes the proof of the 
lemma. □ 



In the light of this we are led to the following conclusion. The three assumptions that 
we made on the functional form of the correction functions are incompatible with the 
conditions for anomaly freedom, unless a, v and a are constants. Therefore, to allow non- 
trivial solution we have to relax one or more of the assumptions which makes the algebra 
much more involved. In the main text of this paper, we organize these calculations by the 
method of counterterms. 
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C Poisson brackets of perturbed variables 



A direct application of the Poisson brackets given by ( 1241) can sometimes be problematic. 
For instance, the Poisson bracket between the two original fields {if, tt}, given by 

{(p + 5ip(x),n + Sir(x)}^ jtt s V M = {<?, tt}^* + {Sip(x), STr(y)} S(p ,sw, 

= l + S(x-y), (154) 

does not agree with the original expression {ip(x), Tf(y)} = 5(x — y). This can be traced to 
the fact that (124")) provides Poisson brackets for the fields (<p, tt, 5ip, 5tt) only if the conditions 
( I2T1) are used in ( l22|) to identify (p and tt with the sole zero modes of inhomogeneous fields. 
The constraints ( 12T1) clearly have a nonzero Poisson bracket {xi,X2}, which makes them 
of the second class. 

According to Dirac [62], second class constraints correspond to non-physical degrees of 
freedom and can be dealt with in the following way. i) One should take linear combinations 
of (all) the constraints, in order to bring as many of them into first class form as possible 
and ii) redefine the Poisson bracket to 

{F, G}s<p,8iT — {F, G}s<p,5n - {F, Xa}s v ,8TvC ab {Xb, G}scp,6ir, (155) 

where 

C ab { Xb ,Xc} = S a c 

so as to remove the variations with respect to the non-physical degrees of freedom. Using 
(I2T]) we obtain 

C n = C 22 = 0, C 21 = -C 12 = (F AiA 2 )-\ 

which implies 

{F,G YsvM = { F,G} SrM - JL (/d 3 ^ / dV^ - (F » o) (156) 

Let us first point out the basic properties of the Dirac bracket (I156p . For the field 
perturbations 

{6<p(x), S7r(y)}*^ s H = 8(x - y) - — . 

Clearly the last term would remove the extra contribution in (11541) yielding the expected 
result 

Thus the Dirac bracket ensures a correct transition from the full theory to the perturbed 
one. By construction, the constraints (I2T!) now commute, {xi>X2}* = 0, and we can 
impose xi an d Xi strongly. Moreover, Dirac brackets between a first order functional and 
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a functional of arbitrary order vanish, which can be seen by inspection using (1156]) . Thus 
for any two functionals, their linear terms do not contribute to the Dirac bracket 

We will now address an issue directly related to closure of the constraints algebra. 
When computing Poisson brackets in the context of perturbation theory, one has a choice 
between two methods: 

1) Calculate the Poisson bracket of the constraints with respect to the full fields and 
expand the resulting expression in orders of perturbations. 

or 

2) Expand the constraints first and then compute their Poisson (Dirac) brackets in 
terms of the expanded fields. 

It is, in general, not guaranteed that the two approaches agree for arbitrary functionals 
which depend on the fields and their first derivatives. However, we have 

Lemma 3 If the fields (p, tt, Sip and 5tt enter the functionals 

F = yd 3 x/(y?,Vy?,7r,V7r), G = Jd 3 xg(ip, V<^, tt, Vtt) (157) 

only as a combination ip = tp + 5ip or tt = tt + 5tt, then the two procedures described above 
yield the same result for {F, G}. 

Proof: We shall show that the second procedure is equivalent to the first one. First of all, 
as linear terms do not contribute to (I156p . we can rewrite the Dirac bracket between two 
expanded constraints as 

{F (o) + F ( 2)) G (o) + G my_^ M = {F (o) + F (D + F W G (o) + G (i) + G my__^ M 

= \ F 1 ^} tp,Tr,5ip,6n' 

According to (I156p . we have 

1 fdFdG _dFdG\ f 3 / 5F 5G 5F 5G 
{ - y Q y-gp-fof ~ dTTd^J + J d x \5{5p>) 5{5tt) ~ 5{5tt) 5{Sip) 

1 f J3 9f f j3 ,dg f j3 8F8G . 



f. 3 SFSG . , 



□ 



6 We are mostly interested in linearized equations of motion, i.e. in the functionals (constraints) we 
should keep terms up to the second order. 
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In the second equality, we have used 



df df 5F S_F = d£_ d _ df 



dip dip' 1 S(5ip) Sip dip d(d a <p) 

and dropped the surface integrals originating from integration of the total divergence terms. 
It is now easy to see that if higher derivative terms were present in the functionals, they 
would have merely led to additional surface terms and would not have affected the final 
conclusion. 

Since linear functionals do not contribute to Dirac brackets, they can be omitted, and 
one can restrict consideration to terms of the zeroth and second order only. Moreover, for 
functionals of an even order, the second term in the Dirac bracket (11561) vanishes, and one 
can simply use the Poisson bracket (f!M|) . 

A somewhat similar consistency issue arises when it comes to equations of motion, 
generated e.g. by a (Hamiltonian) constraint 

H = Jd 3 xh(p, Vp,7t, Vtt). (158) 

There are again two approaches: i) either derive the equations of motion for the original 
fields and then split into the background and (linear) perturbation or ii) expand the con- 
straint and obtain separately equations of motion for the homogeneous and inhomogeneous 
parts of the field. In other words, one needs to compare 

{ip,H} with {<p, H}* and {5p, H}* . (159) 

We start by noting that 

whereas the equation of motion for the background field 
ip = W,H}^ = - / d 3 x— = - / dV 



Vq J dn V J dn 
1 f ^ f dh o dh \ 1 [.» 5h 



~ V J dX {dK da d(d a n) ) ~ V J dX 5ir 

coincides with the background part of the equation of motion (11601) for the total field, 
(5H/8irp \ At the same time, the equation of motion for the perturbation 

5H If,, 5H 



5(5n(y)) 
j d 3 y (^5(x-y)-y^j 



5H /5H\ {1) 
5n(y) \5ti J 
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is nothing else but the perturbed part of the equation (11601) . In fact, one can think of the 
kernel S(x — y) — 1/Vq as cutting off the background part of the function, with which it 
is integrated. It is again pertinent to mention that the linear (as well as the background) 
part of the functional does not contribute to the perturbed equation of motion, that is 

6<p(x) = {6<p(x), H}* SipM = {Scpix), H^} S ip,Sn • (161) 

Note that in the second equality the Poisson bracket is used, not the Dirac bracket. 

To summarize, we have shown that in order to proceed to the perturbation theory, the 
Dirac bracket (11561) in terms of the background and perturbed variables should be used. 
Nevertheless, when dealing with already expanded functionals, containing only even order 
terms, the Dirac bracket reduces to the Poisson bracket ff24|) . 

So far we have considered perturbations of a scalar field. Generalization to tensorial 
fields is rather straightforward for any rank. In particular, we need the canonical pair of 
loop quantum gravity, i.e. the extrinsic curvature and densitized triad whose perturbations 
have Dirac brackets 

l^iw = {F - G} ^- " 4 (/ dSldV Jpfe)) I(MW) - (F " G) ) • 

(162) 

where F and G are arbitrary functionals of K l a and Ef. 

Of interest is also a generalization of the Dirac brackets to the case of local second-class 
constraints. Let us split the triad and extrinsic curvature into the diagonal and traceless 
parts 

K = < + <* E? = e6? + e?, (163) 

such that 

Xl :=tr^ = 0, X2 :=tr< = 0. (164) 

It is easy to see that the pairs (k, e) and (K l a ,ef) are symplectically orthogonal. Indeed, 
the symplectic structure takes the form 

J d 3 x(e5: + e\){k51 + <) = J d 3 x(3^ + 4<). 

However the constraints (11641) are second class under the tentative Poisson bracket 
rn ^ 1 f j3 {5F5G SF5G\ f j3 f SF SG 5F 5G\ /1/MrN 

Specifically, 

{X2(x), Xi{y)}ni,e- = 35(x- y). (166) 
As before, we define the Dirac brackets 

{F, G}* = {F, G}-j d 3 zd 3 z'{F, X a(z)}C ab (z, z'){ Xb (z), G}, (167) 
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where the matrix C ab (x,y) is now space-dependent and satisfies 

J d 3 yC ab (x,y){ Xb (y),Xc(z)} = 5 a c 5(x-z). (168) 

Using the constraints ( 11641) in the equation above, we find that 

C n (x,y) = C 22 (x,y) = (169) 
C l2 (x,y) = -C 21 (x,y) = ±5(x-y). (170) 

Therefore the Dirac bracket reads 

It is easy to see by inspection that the constraints (11641) indeed commute under this 
Dirac bracket, {X11X2}* = 0, so these constraints may be imposed strongly Also, the 
Dirac bracket between the original canonical variables has the correct expression 

{K{x),E){y)Y R ^ £l = 5 b a 5}5(x - y) = {K(x), E b (y)} K ^. 

Earlier on we have seen that one can still use the Poisson bracket rather than the corre- 
sponding Dirac bracket if one removes from the original constraints (that is before splitting 
the canonical variables) all the terms proportional to the second class constraints arising 
because of the splitting. This still holds for local second class constraints. In the case at 
hand, as soon as all the terms containing traces of the extrinsic curvature and the den- 
sitized triad are omitted, the remaining constraints form the correct algebra under the 
Poisson bracket (11651) . 
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